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Abstract 



This thesis consists of three parts. In the first part we review the quantization 
of Yang-Mills theories and perturbative quantum gravity in curved spacetime. 
In the second part we calculate the Feynman propagators of the Faddeev- 
Popov ghosts for Yang-Mills theories and perturbative quantum gravity in the 
covariant gauge. In the third part we investigate the physical equivalence of 
covariant Wightman graviton two-point function with the physical graviton 
two-point function. 

The Feynman propagators of the Faddeev-Popov ghosts for Yang-Mills 
theories and perturbative quantum gravity in the covariant gauge are infrared 
(IR) divergent in de Sitter spacetime. We point out, that if we regularize 
these divergences by introducing a finite mass and take the zero mass limit 
at the end, then the modes responsible for these divergences will not con- 
tribute to loop diagrams in computations of time-ordered products in either 
Yang-Mills theories or perturbative quantum gravity. We thus find effec- 
tive Feynman propagators for ghosts in Yang-Mills theories and perturbative 
quantum gravity by subtracting out these divergent modes. 

It is known that the covariant graviton two-point function in de Sitter 
spacetime is infrared divergent for some choices of gauge parameters. On 
the other hand it is also known that there are no infrared problems for the 
physical graviton two-point function obtained by fixing all gauge degrees 
of freedom, in global coordinates. We show that the covariant Wightman 
graviton two-point function is equivalent to the physical one in the sense 
that they result in the same two-point function of any local gauge-invariant 
quantity. Thus any infrared divergence in the Wightman graviton two-point 
function in de Sitter spacetime can only be an gauge artefact. 
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Chapter 1 
Introduction 



Quantum field theory in de Sitter spacetime has been actively studied re- 
cently due to its relevance to inflationary cosmologies [lj. Furthermore, the 
current observations indicate that our Universe is expanding in an acceler- 
ated rate and may approach de Sitter spacetime asymptotically [2J. The 
original part of this thesis consists of two parts. In the first part we find 
effective Feynman ghost propagators for Yang-Mills theories and perturba- 
tive quantum gravity in de Sitter spacetime. In the second part we find 
the equivalence of the covariant Wightman graviton two-point function with 
arbitrary value of gauge parameters to a infrared finite physical Wightman 
graviton two-point function. 

In order to study higher-order quantum effects for Yang-Mills theories 
or perturbative gravity, one needs to introduce Faddeev-Popov ghosts ex- 
cept in unwieldy gauges such as the axial gauge in Yang-Mills theories [3]. 
In this thesis we shall study the Feynman propagators for the ghosts in 
Yang-Mills theories and perturbative quantum gravity in de Sitter space- 
time. We shall find that ghosts and anti-ghosts for Yang-Mills theories sat- 
isfy the minimally-coupled massless scalar field equation. It is known that 
there is no Feynman propagator for minimally-coupled massless scalar fields 
that respects de Sitter invariance [I], due to infrared (IR) divergence of the 
Feynman propagator for these fields in de Sitter spacetime. We also show 
that the Feynman propagator for ghosts in perturbative quantum gravity is 
also infrared divergent in de Sitter spacetime. Thus if the ghost fields in these 
theories were physical fields, we would need to break de Sitter invariance of 
the vacuum for these fields However we shall see that this problem can 
be circumvented because they are unphysical fields. The interaction between 
the Yang-Mills/gravitational field and the ghosts is such that, if we regu- 
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larize the infrared divergences by introducing a small mass term, the modes 
responsible for the infrared divergences will not contribute in the computa- 
tion of time-ordered products of physical fields. For this reason, we propose 
that one should regularize the infrared divergences of the Feynman propaga- 
tors for the ghosts and anti-ghosts in these theories and then take the limit 
where the regularization is removed. This proposal is equivalent to using 
the effective Feynman propagators obtained by subtracting the regularized 
modes responsible for the infrared divergence in perturbative calculations. 

In fact this effective Feynman propagator for Yang-Mills theories has 
been derived and used in a different context. It was used in calculating 
the covariant graviton propagator two-function in de Sitter spacetime [6] . It 
may be noted that even though the physics there is very different from the 
physics we are considering, the mathematics involved is the same for both 
these cases. 

Infrared divergences in graviton two-point functions have been a matter 
of contention for over two decades [7] . Since linearized gravity has gauge in- 
variance, it is important to determine whether or not these IR divergences are 
gauge artefacts. The graviton two-point function obtained by fixing all the 
gauge degrees of freedom is called the physical graviton two-point function. 
The physical graviton two-point function in spatially flat coordinate system 
was analyzed in Ref. [8]. It was found that this two-point function is IR 
divergent. Infrared divergence of the physical graviton two-point function in 
spatially flat coordinate system was further studied in Ref. [5]. The physical 
graviton two-point function in global coordinate system was obtained in Ref. 
[10] . It was found that this physical graviton two-point function is IR finite. 
The physical graviton two-point function in hyperbolic coordinate system 
was also found to be IR finite [II]. Thus the IR divergence of the physical 
graviton two-point function in spatially flat coordinate system occurred due 
to the coordinate system used and thus was not a physical effect. However we 
will not study the IR divergences of the physical graviton two-point functions 
in this thesis. 

In this thesis we will present a work relevant to the IR divergences in the 
graviton two-point functions that occur for some choices of gauge parameters. 
Allen found that in the covariant gauge the graviton two-point function has 
IR divergences for certain values of the gauge parameters [12] . Thus if a 
gauge is chosen with one of these values of gauge parameters, we shall get 
IR divergences. The graviton two-point function used by Antoniadis and 
Mottola [T3] was IR divergent as they used one of these values of gauge 
parameters. Covariant graviton two-point function with different values of 
gauge parameters was obtained in Ref. [6]. Covariant graviton two-point 
function with arbitrary values of gauge parameters was obtained in Ref. [H] . 
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In this thesis we shall show that the covariant Wightman graviton two-point 
function with any choice of gauge parameters is physically equivalent to the 
physical one obtained in Ref. [TO] . 

For this purpose, we will split the covariant Wightman graviton two- 
point function into vector, tensor and scalar parts [14J. It will be shown 
that apart from a part of tensor part, that also contributes to the physi- 
cal graviton two-point function, all the contributions are pure gauge in the 
following sense. A contribution to a Wightman graviton two-point function 
will be said to be pure gauge if at least in one of its two slots it is a sum 
of tensors obtained by symmetrized derivatives acting on vectors. So a pure 
gauge contribution will be of the form: V a K b 1 a , b ,(x, x') + Vf,-K"^ a / b /(x, x') + 
^ a'Kl bb ,(x,x') +V v Kl ba ,(x,x') for some Kl a , b ,(x,x') and K% a , bl (x, x'). The 
two-point function of a local gauge-invariant tensor field will be the same for 
two Wightman two-point functions if they differ from each other by a pure 
gauge contribution. The scalar part has in fact already been shown to be 
pure gauge in a restricted gauge in Ref. [17]. In this thesis we show that 
this is true for any value of the gauge parameters. We also show that the 
contributions coming from the vector part and a part of the tensor part are 
also pure gauge. Thus the only non-gauge contribution will come from the 
other part of the tensor part. This contribution will be shown to be exactly 
equal to the physical graviton two-point function obtained in Ref. [ID] . 

This physical graviton two-point function obtained in Ref. [ID] suffers 
from no IR divergences. Now if any graviton Wightman two-point function 
in linearized gravity in de Sitter is physically equivalent to that obtained in 
Ref. [ID], as we show in this thesis, then any IR divergences in the covari- 
ant graviton Wightman two-point function in linearized gravity in de Sitter 
spacetime for example as noticed in Ref. [T3], has to be a gauge artefact. 
Thus the IR divergences will not show up in the graviton Wightman two- 
point function of any local gauge-invariant tensor field, e.g., the linearized 
Weyl tensor [15J . 

The rest of the thesis is organised as follows: We shall review Yang-Mills 
theories and perturbative quantum gravity in Chapters 2 and 3, respectively. 
Then after reviewing some basic properties of de Sitter spacetime in Chapter 
4 we shall go on to study ghost fields in de Sitter spacetime. We shall derive 
the explicit expression for the effective ghost propagator for Yang-Mills theo- 
ries in Chapter 5 as a warm up for the derivation of the explicit expression for 
the effective ghost propagator for perturbative quantum gravity in Chapter 
6. Then we shall describe a general method for calculating the Wightman 
two-point function for free field theories in Chapter 7. We shall apply this 
formalism to linearized quantum gravity in Chapter 8. Finally we shall show 
explicitly that the covariant graviton Wightman two-point function with any 
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choice of gauge parameters is physically equivalent to the physical one ob- 
tained in Ref. [10J in Chapter 9. We shall conclude the thesis in Chapter 



We shall do our calculations for Feynman ghost propagators on S A and 
then analytically continue the result to de Sitter spacetime in Chapters 5 and 
6. However we shall do our calculations directly in de Sitter spacetime while 
calculating the Wightman two-function in Chapters 7 and 9. So in Chapter 
5 and 6 propagator means a Feynman propagator obtained by analytical 
continuation from S A and in Chapters 7 and 9 a two-point function means a 
Wightman function calculated directly in de Sitter spacetime. 

The signature we adopt is (— , +, +, +) and we denote a spacetime point 
as x, the spatial component of x as x and its temporal component as t. We 
define the Riemann tensor, Ricci tensor and scalar curvature by 
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[v b v a - v a v 6 ; 



(i.i) 

(1.2) 
(1.3) 



Rbd 

R 




where 



V b A c 




d b A c - T d bc A d , 

^9 da [db9ac + d c g ab - d a g bc ] 



(1.4) 
(1.5) 
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Chapter 2 

Yang-Mills Theories 



In this chapter we shall discuss the quantization of Yang-Mills theories in 
general curved spacetime. After reviewing the quantization of Yang-Mills 
theories we shall also discuss the BRST symmetry for Yang-Mills theories. 



2.1 Basic Formalism 

Yang-Mills theories were originally proposed to explain occurrence of isospin 
[IS] . They have been used to unify the electromagnetism and the weak force 
into a single electro- weak force [19J. Quantum Chromodynamics, which is an 
asymptotically free non-abelian gauge theory, is now thought be the theory 
of strong interaction [20]. Thus except gravity all the forces of nature can 
be described by Yang-Mills theories. Apart from electromagnetism which is 
invariant under £7(1) symmetry all the other gauge theories describing forces 
in nature are invariant under SU(N) symmetry. So we shall first review 
SU(N) Lie groups and some of their properties (see Chapter 15 of Ref. [2T]). 
From now on we shall use the convention that the repeated group indices are 
summed over. An element u of SU(N) can be written as 

u = exp[igA A T A }. (2.1) 

Here T A are N x N traceless Hermitian matrices, which satisfy 

Tr(T A T B ) = U AB , (2.2) 

The matrices T A are the generators of SU(N) Lie algebra, which is given by 

[T A ,T B ]=if% B T c , (2.3) 



9 



where f% B are called the structure constants. Strictly speaking T A form a 
basis of the defining representation of the SU(N) Lie algebra. 

These structure constants are antisymmetric in their lower indices: as 

[T A ,T B ) = -[T B ,T A ], (2.4) 

we have 

f A B = ~Iba- (2-5) 
They also satisfy the Jacobi identity: as 

[T A , [T B , T c }} + [T B , [T c , T A }\ + [T C , [T A , T B )} = 0, (2.6) 

we have 

IadIbc + fsofcA + fcofAB = 0- (2-7) 

We can now define matrices r A which form a basis for the adjoint represen- 
tation of this Lie algebra: as 

{ta) B c = iflc- (2-8) 

Now from the Jacobi identity given in Eq. (12. 7p . r A satisfy the Lie algebra 
given in Eq. (12. 3p . 

[TA,T B }=lf% B T C . (2.9) 

Having reviewed some properties of SU(N) Lie algebra, we shall briefly 
review SU(N) Yang-Mills theories (see Chapter 9 of Ref. [22J). We start 
from a matter field 4>(x) which forms a N- dimensional vector in the defining 
representation space of the Lie algebra. So it transforms under a spacetime 
dependent SU(N) transformation u(x) as follows: 

(f) u (x) = u(x)(f)(x). (2.10) 

We now also want a derivative of <j)(x) that transforms like <p{x). However 
V a (fi(x) does not do so. So we define a derivative D a <p{x) called the covariant 
derivative which transforms like <j)(x). It is defined to be 

D a <f>=\y a -igA a ](f>, (2.11) 

where A a is defined to be a matrix-valued gauge field. It can be given by 

A a = A A a T A . (2.12) 

Now we expect the covariant derivative to transform like <fi{x). So if A u a is 
the transformed gauge field, we have 

[V - igAZ]u<f> = uD a <f>. (2.13) 
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So we have 

\y a -igA«]u(i>-uD a <j> = [{V a u)-igA u a u + iguA a ] ( i ) = 0. (2.14) 
Now we have 

igA u a u = {V a u) + iguA a . (2.15) 
Multiplying by —ig~ x u~ x from the right-hand side, we get 

A u a = uAaU' 1 - ig-\V a u)u-\ (2.16) 

Now we expand u and u^ 1 as 

u = l + tgA A T A + 0(A 2 ), (2.17) 
u- 1 = l-igA A T A + 0{A 2 ). (2.18) 

The gauge field A™ can now be written, to first order in A A : as 

Al = [l + igA A T A ]A a [l-igA B T B ]-ig-\V a igA A T A ) 
= A a + igA A [T A ,T B ]A B a +V a A A T A 
= A a -gfi B T c A A A B a +V a A A T A 

= \A A + V a A A + gf A c A B a A c ]T A . (2.19) 
So if the infinitesimal transformation of A a is given by 

5 A A a = T A 5 A A A , (2.20) 

then we have 

5 A A A = V a A A + gf A c A B A c . (2.21) 



Now we define F A h as follows: 



This is the gauge transformation of the gauge fields. 

[D a ,D b ]<j ) = tg~ 1 F A <j ) , (2.22) 
F ab = F A T A . (2.23) 



where 
Thus we get 



F ab <P = g ^gVaAb- gV b A a -ig 2 [A a ,A b ]](f) 

= T A [V a A A -V b A A + gf A c A B Af]<f>. (2.24) 



So we get 



F A = V a A A - V b A A + gji c A B a A c b . (2.25) 
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Now by construction [D a , D b ](j) transform like <fi: 

[D a ,D 6 ]0->u[A,,Z> 6 ]0. (2-26) 
This implies that transforms as follows: 

F ab = uF ab u-\ (2.27) 
This is because it follows from Eq. (I2.26P 

F ab <p -»■ uF ah u- x u<\) = uF ah <\). (2.28) 
The classical action for the gauge field is given by 

S = I d 4 x^g~C c , (2.29) 



where C c is the classical Lagrangian, which is given by 

£ c = - l -Tr(F ab F ab ). (2.30) 

In this thesis we define the Lagrangian to be a scalar rather than a scalar 
density that would include the measure y/—g. It can be seen to be invariant 
under u transformations, 

Tr(F ab F ab ) Tr(uF ab u- l uF ab u- 1 ) 

= Tr(F ab F ab ). (2.31) 

This classical Lagrangian can also be written as 

1 



because 



C c = — A F*Ff, (2.32) 



C c = ~Tr(F£F abB T A T B ) 



F*F abB 5 AB 



4 

here we have used Eq. (12.21) . 



(2-33) 
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2.2 Gauge Fixing 



In this section we discuss gauge-fixing for Yang-Mills theories. Fixing a gauge 
is essential in Yang-Mills theories before they can be quantized. To see this 
we first shall briefly review the quantization of the scalar field theory. If we 
consider a massive free scalar field on a (3 + 1) dimensional spacetime, then 
the Lagrangian of this scalar field is given by 

£ = ~[-V a 0V a 0-m 2 2 ]- (2.34) 

We can canonically quantize this theory by first defining conjugate momen- 
tum 7r(x) as 

vr = V^j^f- = -V=gg 0a v a <P, (2.35) 

then by imposing the following canonical commutation relations: 

[0(x,t),7T(x',t)] = Z5(X,X'), 

[<j>(x,t),<f>(x',t)) = 0, 

[7r(x,t),7r(x',t)] = 0, (2.36) 
where <5(x, x') is defined as 

J d 3 x5(x,x')/(x) = /(x'), (2.37) 

for any compactly supported smooth function /(x). Note here that (i 3 x = 

doc diOc di^jc > 

To apply the above procedure to Yang-Mills theories, we would have to 
define 

dC 

Tr a A = J=g~- j. (2.38) 

However, the conjugate momentum 7r^ corresponding to A A is constrained 
to vanish: 

So we can not quantize Yang-Mills theories simply by imposing the canonical 
commutation relations, that are the straightforward generalization of those 
used in the scalar field theory. In fact the constraint tt° a = is closely related 
to the fact that the initial Lagrangian given in Eq. (12.321) is invariant under 
gauge transformations [23J, 

6 A A A = V a A A + gf£ c A B a K c . (2.40) 
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However as is well known, we can overcome this problem by fixing the gauge 
and this can be achieved by adding a gauge-fixing term and a ghost term to 
the original classical Lagrangian |24j . 

Let us explain the general procedure for obtaining a Lagrangian which 
can be quantized following the canonical quantization procedure [25]. We 
start by choosing a gauge-fixing function, GL^. In the classical theory one 
imposes the condition GL4] A = 0. In the quantum theory one adds a gauge- 
fixing term and a ghost term to the classical Lagrangian. The gauge-fixing 
term is obtained by first squaring this gauge-fixing function GL4] A in such a 
way that it forms a scalar in the representation space of the Lie algebra and 
then multiplying it by —l/(2a). To find the ghost Lagrangian we first take 
the gauge transformation of the gauge-fixing function and then change 

the A A to ghost fields c A . After that we contract any free index left with 
anti-ghosts c A to form a scalar quantity in the representation space of the 
Lie algebra and then multiply it by i. The scalar thus obtained is called the 
Faddeev-Popov ghost Lagrangian [24J. Both the ghost fields and anti-ghost 
fields obey Fermi-Dirac statistics. The Faddeev-Popov ghost Lagrangian is 
also needed for ensuring the unitarity of the theory (see Chapter 16 of Ref. 
[21 j ) . In fact Faddeev-Popov ghosts were first conjectured by Feynman using 
unitarity [26] and later derived by Faddeev and Popov by using path integral 

In the classical theory, in Lorentz gauge the gauge-fixing condition is 
given by 

G[A] A = S7 a A Aa = 0. (2.41) 

In the quantum theory the gauge-fixing term C g corresponding to this gauge- 
fixing condition is obtained by first squaring this gauge-fixing condition in 
such a way that it forms a scalar quantity and then multiplying it by —1/ (2a). 

C g = -^[V a A\V b A Ab }. (2.42) 

To obtain the ghost Lagrangian we first take the gauge transformation of the 
gauge-fixing function: 

5 A G[A} A = V a [V a A A + gf^cK^l (2-43) 

and then change A A to c A , 

5 {C) G[A] A = V a [V a c A + gf A c A B c c ). (2.44) 

Finally we contract with anti-ghosts c A , in such a way that it forms a scalar 
quantity, 

c A 5 {c) G[A] A = c A V a [V a c A + gf A c A B a c c ). (2.45) 
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After dropping total derivative terms, and multiplying by i we obtain the 
ghost Lagrangian, 

C gh = -iV a c A D a c A , (2.46) 

where D a c A is the covariant derivative for the fields in adjoint representation 
with generators t a , 

D a c A = V a c A + gf A c A B c c . (2.47) 

The total Lagrangian is given by the sum of the classical Lagrangian, the 
gauge-fixing term and the ghost term, 

C = C c + C g + C gh . (2.48) 



2.3 BRST Symmetry for Yang-Mills Theo- 
ries 

The original classical Lagrangian is invariant under gauge transformations. 
After fixing the gauge this invariance is broken. However the total La- 
grangian, which is obtained as a sum of the original classical Lagrangian, 
the gauge-fixing term and the ghost term is now invariant under a transfor- 
mation called the BRST transformation [27] . To see this we first rewrite the 
gauge-fixing term in terms of an auxiliary field B A , as: 



4, 



B A V a A Aa + ^B A B A 



(2.49) 



It is possible to recover the original gauge-fixing term from this using the 
field equation for B A . Now the BRST transformation for Yang-Mills theory 
is given by 

s = es, (2.50) 
where e is an anti-commuting complex number and 



sA A = i[V a c A + gf A c A B a c c ] 
sB A = 0, 
sc A = B / 



sc A = -jf A c c B c c . (2.51) 

As the BRST transformation is fermionic in the sense that it depends on 
a fermionic parameter e, it is important to fix a convention as to how the 
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transformation s acts on the fields. We let them act from the left on the 
fields, for example 

s(c A c A ) = s(c A )c A - c A sc A . (2.52) 
Now these transformations can be seen to satisfy 

(s) 2 = 0. (2.53) 

The action of (s) 2 on B A and c A trivially vanishes. 

s(sB A ) = 0, (2.54) 

and 

s(sc A ) = sB A = 0. (2.55) 

Its action on c A is given by 

(sfc A = -js[f A c c B c c } 

= 0. (2.56) 

We have used the Jacobi identity given by Eq. (12.71) and the fact that c A are 
anti-commutating in the last equality. Similarly we have 

{sfA A = is[V a c A + gf A c A B a c c ] 

= {f^a{c B c c )-gf A c iy a c B )c c 

~9 2 fBcfEF^a cF ° C + ~YfBcfEFA B C E C F 

= 0, (2.57) 
where we have used the Jacobi identity given by Eq. (12.71) and the fact that 

f£cVa(c B C C ) = fi C ^aC B )c C + f A c C B V a C C 
= f£c(VaC B )c C -f A B C B V a C C 

= f£c(V a c B )c c -f A c c c V a c B 

= 2f A c (V a c B )c c . (2.58) 

The original classical Lagrangian is invariant under the BRST transforma- 
tion, as the BRST transformation for the original classical Lagrangian is just 
the gauge transformations with A A replaced by c A : 

sC c = 0. (2.59) 
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We also note that 



CA 

B a 



y a A Aa + -B A 



V a A Aa + -B J 
a 



B A V a A Aa + -B A B A 



+ ic A V a D a c A 
- iV a c A D a c A 



C + C 



ghi 



(2.60) 



up to a total divergence term. The invariance of the sum of the gauge-fixing 
term and the ghost term can now be deduced from the fact that their sum 
is written as a total BRST variation. Thus the BRST variation of the sum 
of the gauge-fixing term and the ghost term vanishes due to Eq. (12.53p . as: 



s[C g + C 



gh\ 



CA 



V a A Aa + -B J 



a 



0. 



(2.61) 



However, as the original classical Lagrangian is also invariant under the 
BRST transformation, so the total Lagrangian will also be invariant under 
the BRST transformation: 



S[C C + Cg + C 



gh\ 



0. 



(2.62) 



The invariance of Yang-Mills theory under BRST symmetry is crucial in 
selecting the physical states. 



2.4 Physical States 

In Abelian gauge theory in flat spacetime, we get negative norm states in 
the Fock space. For example, for free Abelian gauge theory in the Feynman 
gauge the creation and annihilation operators, a a (k) and a\{k') satisfy, 

[a a (k),al(k')}=tg ab 5(k,k'), (2.63) 

so we have 

[a (k),al(k')] = -id(k,k'), (2.64) 

It is possible to restrict our Fock space to positive norm states in Abelian 
gauge theory by using Gupta-Bleuler formalism [28]. In this formalism phys- 
ical states \(f>p) are defined by requiring that 

d a A^ a (x)\<P P ) =0, (2.65) 
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where d a A^ a (x) is the positive frequency part of d a A a {x). It can be shown 
that 

Dd a A a (x) = 0, (2.66) 

even if the field A a interacts with charged matter. Thus, this method works 
for Abelian gauge theory in Minkowski spacetime, even if interactions are in- 
cluded. However it does not work for non-abelian gauge theories in Minkowski 
spacetime as d a A a (x) does not satisfy the free field equation for non-abelian 
gauge theories and hence one cannot define its positive frequency part. It 
also fails in case of Abelian gauge theory in general curved spacetime as it is 
not always possible to define the positive frequency part of V a A a (x) uniquely 
in general curved spacetime. 

However it is possible to remove these negative norm states in case of 
non-abelian gauge theory in Minkowski spacetime by the Kugo-Ojima crite- 
rion [25J. Kugo-Ojima criterion states that the physical states \<f> p ) must be 
annihilated by Q [25], 

Q\<Pp) = 0, (2.67) 

where Q is the BRST charge, which is the Noether charge corresponding 
to invariance of the total Lagrangian under the BRST transformation. The 
Kugo-Ojima criterion divides the Fock space into two parts, the physical part 
which is annihilated by Q and the unphysical part which is not annihilated 
by Q. The physical part in turn has a trivial part which is composed of 
those states which are obtained by the action of Q on states that are not 
annihilated by Q, \<f) t ) = Q\cj) up ), where Q\4> U p) 7^ 0. These trivial physical 
states are orthogonal to all physical states, 

{(j> p \(j) t ) = (0 P |Q|<M = 0. (2.68) 

Thus in Minkowski spacetime the only relevant physical states for non- 
abelian gauge theory are those that are annihilated by Q and are not obtained 
by the action of Q on any other state. 

In curved spacetime it is expected that BRST symmetry will play an 
important role in defining the physical states [2H]. However, in this thesis 
we shall show that the Feynman propagator for ghost fields in de Sitter 
spacetime suffers from IR divergences. Due to this problem the definition of 
BRST symmetry for Yang-Mills theories in de Sitter spacetime might become 
non-trivial. However if BRST symmetry is defined for Yang- Mills theories in 
de Sitter spacetime then the Kugo-Ojima criterion can be used to obtain the 
physical subspace, in analogy to what is done in case of Minkowski spacetime. 
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Chapter 3 

Perturbative Quantum Gravity 



In this chapter we shall discuss the quantization of perturbative quantum 
gravity in general curved spacetime. We shall also discuss the BRST sym- 
metry for perturbative quantum gravity in this chapter. 



3.1 Basic Formalism 

According to General theory of Relativity, gravity is caused by the curvature 
of spacetime [30]. The Lagrangian for gravity with a cosmological constant 
A, is given by 

£ c = >£iJ[J2</>-2A]. (3.1) 
lo7rG 

( f) 

We have denoted the full metric as g^ h to distinguish it from the fixed back- 
ground metric g ab . The scalar curvature corresponding to this full metric is 
denoted by R^>. We adopt units such that 

16ttG = 1. (3.2) 

The Lagrangian given by Eq. (13.1 j) is invariant under the following infinites- 
imal transformations originating from its general coordinate invariance: 

= £a9 ( £, (3-3) 
where the Lie derivative £\t ab for any tensor t ab is given by 

£ A t ab = A c V c t ab + t ac V b A c + t bc V a A c . (3.4) 
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In perturbative gravity one writes the full metric in terms of a fixed 
background metric and small perturbations around it. We also denote the 
small perturbation around the fixed background metric as h a b- So we can 
now write, 

Gab = 9ab + hab. (3.5) 

This small perturbation is regarded as a field that is to be quantized. The 

covariant derivative along with the raising and lowering of indices will be 

with respect to the background metric. 

For perturbative quantum gravity, we first write the Lagrangian in terms 
f f 1 ( f) 

of the full metric g ab and then expand g ab in terms of h a b and g a b- So we 

get the gravitational Lagrangian for perturbative quantum gravity in terms 

of h a b. It will in general contain infinitely many terms. Each term containing 

two derivatives of h a b comes from the Ricci scalar Rv\ and there are only 

finitely many terms of this kind. However, because the inverse of the metric 

contains infinitely many terms with integer powers of h a b, we have infinitely 

many terms in the Lagrangian for perturbative quantum gravity. 

Now as g a b is fixed, the transformation of g^ b will be attributed to h a b- 

Thus, the transformation of h a b is now given by 



S\h a b = £\g ab 

= £k9ab + £\h a b 

= V a A b + V b A a + £\h ab . (3.6) 

In the last line we have used V c g a b = 0. To first order in A a , the Lagrangian 
for perturbative quantum gravity will be invariant to all orders in h a b, under 
the following transformation, 

hKb = V a A fe + V b A a + £ A h ab , (3.7) 

where the Lie derivative £hh a b is given by 

£ A h ab = A c V c h ab + h ac V b A c + h bc V a A c . (3.8) 

Thus gravity is somewhat analogous to Yang-Mills theories. Just as in the 
Yang-Mills theory the Lagrangian was invariant under gauge transformations, 
here the Lagrangian is invariant under general coordinate transformations 
given by Eq. (13. 7p . Due to this analogy one can deal with gravity in a 
similar way to what was done for Yang-Mills theories. 
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3.2 Quantization 



In the previous section we found that the Lagrangian for perturbative quan- 
tum gravity is invariant under the following transformations: 

5 A h ab = V a A b + V b A a + £ A h ab . (3.9) 

In the previous chapter we found that the Lagrangian of the Yang-Mills 
theory was invariant under a gauge transformation which led to constraints 
and we were not able to quantize this theory without fixing a gauge. Now 
as the Lagrangian for perturbative quantum gravity is invariant under the 
above mentioned transformations, we expect that constraints will exist for 
perturbative quantum gravity also. In fact it is known that just like in the 
Yang-Mills case the invariance of the Lagrangian for perturbative quantum 
gravity leads to the existence of constraints [3 1 J . We are thus not able to 
quantize it using canonical commutation relations without fixing the gauge. 
Thus we add a gauge- fixing term and a ghost term to the classical Lagrangian 
of perturbative quantum gravity, just like what we did in case of Yang-Mills 
theories (see Chapter 2 of Ref. [32]). We shall explain this procedure below. 
In classical theory we impose the gauge-fixing condition, 

G[h] a = (V b h ab - W a h) = 0. (3.10) 

where, 

k^l. (3.11) 

For k — 1, the conjugate momentum for /ioo still vanishes, so we take k ^ 1. 
For this reason sometimes k is written as 1+/3 -1 , where /3 is an arbitrary finite 
constant [14] . In quantum theory we add a gauge-fixing term and a ghost 
term to the classical Lagrangian. The gauge-fixing term corresponding to this 
gauge-fixing condition is obtained by first squaring this gauge-fixing function 
G[h] a , in such a way that it form a scalar quantity and then multiplying it 
by-l/(2«): 

Cg = -^Whab ~ kV a h] [V c h aC - Wh] . (3.12) 

Now to get the ghost action, we follow the same procedure we followed for 
getting the ghost action in the case of Yang-Mills theory. First we take the 
gauge transformation of the gauge-fixing function: 

5 A G[h] b = V a [V a A b + V b A a - 2kg ab V c A c 

+£ A h ab - kg ab g cd £ k h cd \. (3.13) 
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Next we introduce the Faddeev- Popov ghost field c a , which is a fermionic 
vector field, and define £Ji a b to be 

£cKb = c c W c h ab + h ac V b c c + h bc V a c c . (3.14) 

Then we change A a to c a in Eq. (13.131) 

6 (c) G[h] b = V a [V a c b + V b c a - 2kg ab V c c c 

+£ c h ab - kg ab g cd £ c h cd \. (3.15) 

Now we contract the free index with anti- ghost field c , which is also a 
fermionic vector field: 

c b S (c) G[h} b = c b V a [V a c b + V b c a - 2kg ab V c c c 

+£ c h ab - kg ab g cd £ c h cd \. (3.16) 

Finally we multiply it by i and drop total derivatives to get the ghost La- 
grangian, 

C gh = -iV a c b [V a c b + V b c a - 2kg ab V c c c + 
c c V c h ab + h ac V b c c + h bc V a c c - 
kg ab 9 cd [c e V e h cd + h ce V d c e + h de V c c e ]]. (3.17) 

The total Lagrangian is given by the sum of the classical Lagrangian, the 
gauge-fixing term and the ghost term: 

C = C c + C g + C gh . (3.18) 



3.3 BRST Symmetry for Perturbative Quan- 
tum Gravity 

The original classical Lagrangian is invariant under gauge transformations. 
However after fixing the gauge this invariance is broken. However the total 
Lagrangian for the perturbative quantum gravity, which is obtained as a sum 
of the original classical Lagrangian, the gauge-fixing term and the ghost term 
is also invariant under BRST symmetry [33]. To see this we first rewrite the 
gauge-fixing term in terms of an auxiliary field B a , as: 



C g = B a 



V b h ab - kV a h + |s a 



(3.19) 
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The original form of the gauge- fixing term can be obtained from this one using 
the field equation for B a . Now the BRST transformation for perturbative 
quantum gravity is given by 

s = es, (3.20) 
where again e is an anti-commuting complex number and 

sh ab = i[V a c b + V b c a + c c V c h ab 

+h ac V b c c + h bc V a c c ], 
sB a = 0, 



sc a = -ic b V b c a . 
This transformation also satisfies (see Chapter 5 of Ref. [34|) 

(s? = 0. 



(3.21) 



(3.22) 



The original classical Lagrangian is again invariant under the BRST trans- 
formation, as the BRST transformation for the original classical Lagrangian 
is just the gauge transformation with A a replaced by c a : 



sC c = 0. 

We also note that up to a total divergence, 



(3.23) 



s 



C" 



a: 



V b h ab - kV a h + -B, 



a 



V b h ab - kV a h + -B„ 



+ 



irV b [V a c b + V b c a - 2kg ab V c c c 
+c c V c h ab + h ac V b c c + h bc V a c c - 
kg a b9 cd [c e V e h cd + h ce V d c e + h de V c c e ]} 

£-9 + ^gh- 



(3.24) 



The invariance of the gauge-fixing term and the ghost term can be now 
deduced from the fact that their sum is written as a total BRST variation 
and that the BRST variation of a total BRST variation vanishes: 

,2 



s[C g + C 



ghl 



c" 



W b h ab - kW a h + y B a 



0. 



(3.25) 



However as the original classical Lagrangian was also invariant under the 
BRST transformation, so the total Lagrangian will also be invariant under 
the BRST transformation: 



S[C C + Cg + C 



gh\ 



0. 



(3.26) 
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The invariance of perturbative quantum gravity under BRST symmetry is 
also crucial in selecting the physical states of the theory It is possible to 
define physical states in perturbative quantum gravity, at least formally using 
Kugo-Ojima criterion (see Chapter 5 of Ref. [S]). Thus, the state \4> p ) is 
said to be a physical states if 



The only relevant physical states in perturbative quantum gravity again are 
those states which are annihilated by Q and are not obtained by the action 
of Q on any other state. 

In this thesis we will show that the Feynman propagator for ghost fields 
in perturbative quantum gravity in de Sitter spacetime also suffers from IR 
divergences. Thus just like the Yang-Mills case, due to IR divergences the 
definition of BRST symmetry for perturbative quantum gravity in de Sitter 
spacetime might become non-trivial. 




(3.27) 
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Chapter 4 

de Sitter Spacetime 



In this chapter, we shall review some basic properties of de Sitter spacetime 
(see Chapter 5 of Ref. [30]). We shall start by reviewing properties of dif- 
ferent tensors in de Sitter spacetime and then present some basic coordinate 
systems used in de Sitter spacetime. 



4.1 Properties of de Sitter Spacetime 

The vacuum Einstein equations with a cosmological constant A are given by 
(see chapter 3 of Ref. [30] ) . 

G ah = -Ag ab , (4.1) 
where G ab is the Einstein tensor which is given by 

Gab = Rab — -Rg a b- (4.2) 

Constant curvature spacetimes are solutions to Eq. (14.21) . and are character- 
ized by the condition (see Chapter 5 of Ref. [30J) 

R a bcd = ^R[5 a c g bd -5 a d g bc ], (4.3) 

where R is a constant. The spacetimes with R = 0, R < and R > are 
called the Minkowski spacetime, the anti-de Sitter spacetime and de Sitter 
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spacetime, respectively. Now for the constant curvature spacetimes, we have 



Rbd — R a bad 

= ^R[K9m - s a d9ba ] 

= \Rg db . (4.4) 

So the Einstein tensor G ab for the constant curvature spacetimes is given by 

Gab = Rab ~ ~Rg ] ab 

= -^3abR — -jRQab 

= -\g ab R. (4.5) 

Now if we compare Eq. (14 .2p with Eq. (14.51) . we see that we can view Eq. 
( 14. 5p as a solution of the vacuum Einstein field equations with a cosmological 
constant A given by 

A = -R. (4.6) 
4 

As de Sitter spacetime is defined to be a spacetime of constant positive 
curvature, so for de Sitter spacetime the cosmological constant is positive. 
This cosmological constant is related to a constant called the Hubble constant 
H, as follows: 

A = 3H 2 . (4.7) 

It is called the Hubble constant as it was used by Hubble for the measure 
of the expansion of the universe [SB] ■ So for de Sitter spacetime, in terms of 
the Hubble constant if, we have 



R = 12H 2 , (4.8) 

R ab = 3H 2 g ab , (4.9) 

R a bcd = H 2 [5 a c g bd -5 a d g bc ], (4.10) 

G ab = -3H 2 g ab . (4.11) 



4.2 de Sitter Spacetime Metric 

de Sitter spacetime, which is defined to be a spacetime of constant positive 
curvature, has the topology R x S 3 and can be viewed as a hyperboloid 
in five dimensional Minkowski spacetime (see Chapter 5 of Ref. |30j). If 
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the Cartesian coordinates in the five dimensional Minkowski spacetime are 
X, Y, Z, W, T and the metric is given by 

ds 2 = -dT 2 + dW 2 + dX 2 + dY 2 + dZ 2 , (4.12) 

with — oo < W, X, Y, Z, T < oo then de Sitter spacetime is the hypersurface 
given by the following equation, 

X 2 + Y 2 + Z 2 + W 2 -T 2 = a 2 . (4.13) 

Here a is called the radius of de Sitter spacetime and is related to the Hubble 
constant by a = H^ 1 , as we shall see. De Sitter spacetime inherits the 
five dimensional Lorentz invariance of five dimensional Minkowski spacetime, 
which is 50(4, 1) [37]. This is called de Sitter group and is defined to be the 
group of linear transformations in five dimensions which preserve Eq. f)4. 13|) . 
A convenient parametrisation of this hypersurface is given by 

T = asinh(a _1 i), 

W = acosh(a _1 t) cos-0, 

X = acosh(a _1 t) sin ■?/> sin 6* cos 0, 

Y = a cosh(a _1 t) sin ■?/> sin # sin 0, 

Z = acosh(a _1 t) sin -0 cos ^. (4-14) 

with < 0, 6 < 7r, < <p < 2tt and — oo < t < oo. Then substituting Eq. 
( I4.14p in the Minkowski metric given by Eq. (I4.12p . we get 

ds 2 = -dt 2 + a 2 cosh 2 a -H[dip 2 + sin 2 ip(d6 2 + sin 2 6d<f) 2 )}. (4.15) 

The singularities at ip = 0, tt and 9 = 0, it are those singularities that occur 
in polar coordinates. Apart from them this coordinate system is regular in 
the whole of de Sitter spacetime and covers all of de Sitter spacetime. If we 
rescale by letting t — > at, then the metric given in Eq. (14 . 1 5 j) becomes 

ds 2 = a 2 [-dt 2 + cosh 2 t[dip 2 + sin 2 iP(d9 2 + sin 2 0# 2 )]]. (4.16) 

If we further perform the transformation, 

X=\-*, (4-17) 
then this rescaled de Sitter metric given in Eq. (I4.16P becomes 

ds 2 = a 2 [d X 2 + sm X 2 [d^ 2 + sin 2 ^[d9 2 + sin 2 0# 2 )]]. (4.18) 
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This is in fact the metric on a four dimensional sphere denoted by S* 4 , of 
radius a. It may be noted that for \ to be real in Eq. ( 14.1 7p . time t has to 
be imaginary. So we are using imaginary time in Eq. ( 14.181) . 
We can also write the de Sitter metric as follows: 



ds 2 = -dt 2 + exp(2a~H)(dx 2 + dy 2 + dz 2 ), (4.19) 



where 



(W + T) 
t = a log . 

a 

X 

x = a— — — , 
W + T 

Y 

v = a- 



W + T' 
Z 

z = a— — -. 4.20 
W + T K ' 

Here — oo < x, y, z < oo and < t < oo. However these coordinates only 
cover half of de Sitter spacetime as t is not defined for W + T < 0. 
Now if the FRW metric is given by 

ds 2 = -dt 2 + a 2 (t)dn 2 , (4.21) 

where dfl 2 is the spatial part of the metric. It can be shown that the Hubble 
constant given in Eq. (14 .7p . is related to a by (see Chapter 5 of Ref . [35] ) 

H = -. (4.22) 
a 

So for Eq. (14. 19p . the Hubble constant is given by 

H = -. (4.23) 

a 

Unless specified otherwise, from now on we shall choose units such that 
H 2 = 1, to simplify calculations. 
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Chapter 5 



Ghosts for Yang-Mills Theories 
in de Sitter Spacetime 



In this chapter we shall examine the ghost propagators for Yang-Mills theories 
in de Sitter spacetime. We shall use the fact the Feynman propagators in 
the Euclidean vacuum [I] in de Sitter spacetime can be obtained from the 
corresponding Green's function on the S A by analytical continuation |38j . 

5.1 Ghosts for Yang-Mills Theories 

We have already presented the ghost Lagrangian for Yang-Mills theories in 
general spacetime in Eq. ( I2.46p . So we can start from the following ghost 
Lagrangian, 

C gh = -iV a c A D a c A . (5.1) 
Now as the covariant derivative D a c A is given by 

D a c A = V a c A + gf£ c A*c c , (5.2) 

we have 

C gh = -i[V a c A V a c A + gf A c V a c A A B a c c ). (5.3) 

This can now be written as a free Lagrangian C^ e and an interaction part 
Oghi where the free Lagrangian is given by 

C f r = -iV a CAV a c A . (5.4) 
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The equation of motion obtained from this free Lagrangian for ghosts and 
anti-ghosts are 

Uc A = 0, (5.5) 
Uc A = 0. (5.6) 

These are the equations of motion for the minimally-coupled massless scalar 
fields. 

As is well known, there is ambiguity in the selection of positive frequency 
solutions and thus an ambiguity in selecting the vacuum state in curved 
spacetime [39]. However, for de Sitter spacetime there is a de Sitter invari- 
ant vacuum state called the Euclidean vacuum [3j in which the Feynman 
propagator can be obtained from the Green's function on the four sphere by 
analytical continuation [SB]- We shall use the Euclidean vacuum state for 
calculating the Feynman ghost propagator for Yang-Mills theories. 

Let |0) be the Euclidean vacuum state. Then the Feynman propagator 
for the free ghost fields would be given by 

(0\T[c A (x)c B (x')]\0) = i5 AB D (x,x r ), (5.7) 

where D (x, x') would satisfy 

nD Q {x,x') = -5\x,x'). (5.8) 

However, it is well known that there is no Feynman propagator for minimally- 
coupled massless scalar fields that respects de Sitter invariance [3]. We can 
see that there is no Euclidean vacuum for these fields as follows. To find 
the Feynman propagator in the Euclidean vacuum, we first find the Green's 
function on S A . Now a complete set of basis for any scalar function on S 4 is 
the scalar spherical harmonics Y La which satisfy [40J, 

- uY La = L(L + 3)Y La J (5.9) 

where L = 0, 1,2,3,4- ■ ■ and a represents all the other labels. The scalar 
spherical harmonics are normalized as: 

J d 4 x^Y La Y* L ' a ' = 5 LL '5 aa '. (5.10) 

We define 8(x,x') as: 

J v^jd 4 x'f(x')5 A (x, x') = f{x), (5.11) 
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for any function f(x'). We can mode expand 6(x, xf) as: 

oo 

5 4 (x, x') = E Y L °(x)Y* La (x'), (5.12) 

L=0 a 

because if we mode expand any function f(x') on S A as: 

oo 

/M-EE^'^'M. (5-13) 

L'=0 cr' 

where Xva 1 are constants, then use Eq. (I5.10p . we find 

oo „ 

EE / V9d*x'f(x')Y L °(x)Y* L °(x>) 

L=0 cr •* 

oo oo „ 

= EE EE / V9d 4 x>\ L , a ,Y LW (x')Y L °(x)Y* L °(x') 

L=0 a L'=0 cr' J 

oo 

= LEa^co = /(*)• ( 5 - 14 ) 



L=0 cr 



If Eq. f)5.8p were to be satisfied, the Green's function D (x, x') would have 
to be decomposed into spherical harmonic modes as follows: 

oo 

D (x, x') = E k L Y L(T (x)Y* La (x'). (5.15) 

L=0 cr 

where is a constant. Now substituting Eqs. (I5.12p and (15.151) into Eq. 

we get, 

oo 

□ E k L Y La (x)Y* La (x') = - E E Y La (x)Y* La (x'). (5.16) 

L=0 cr L cr 

From Eqs. (15.91) and (I5.16p . we get 

fci = Z(zW (5 - 17) 

So formally we can write 
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However, D (x,x') is actually not well defined because L(L + 3) = for 
L = 0. If we regulate this propagator by adding a small mass m 2 to Eq. 
f!5.8p . then we get 

[□ — m 2 ]D m 2(x, x') = — 8(x, x). (5.19) 

Then following what we did for the minimally-coupled massless scalar field, 
we get 

V m2 {x,x)- 2^2^ L(L + 3) + m2 " ^ 2Uj 

L=0 a- v ' 

This propagator given by Eq. (I5.20p diverges in the zero mass limit. Such 
divergences are called infrared divergences. As the volume of the unit S n is 
given by 

27T n+1 

so the volume of the unit S 4 is 

(5.22) 

Now from from Eqs. ( 15. lUj) and (I5.22p . the L = mode will be given by 



The L = mode, which is a constant mode, is the cause of infrared divergence 
as its contribution diverges in the zero mass limit. Thus by substituting Eq 
(E22D into Eq. (1QP|) . we find 

L=l a K ' 

This propagator obviously diverges in the zero mass limit. 



5.2 Effective Propagator 

We saw in the previous section that the propagator for the ghost fields suffers 
infrared divergence in de Sitter spacetime as the free part of the ghost La- 
grangian satisfies the minimally-coupled massless scalar field equation, which 
is known to be IR divergent in de Sitter spacetime. If the ghost fields were 
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physical fields, we would need to break the de Sitter invariance of the vac- 
uum for these fields [5]. However we shall see that this problem can be 
circumvented because they are unphysical fields and thus do not exist in 
the final/initial states of any scattering calculations and they only appear in 
internal loops in the Feynman diagrams. This can be seen as follows. 

First we note that if we consider the propagator which is regulated by 
addition of a small mass, then we observe that the infrared divergence are 
caused by the L = mode, which is a constant mode, in the zero mass limit. 
However this constant mode does not contribute to the loops in the Feynman 
diagrams. This is because the interaction part of the ghost Lagrangian is 
given by 

£3? = -igf^ a c A .A B a c c , (5.25) 
so the anti-ghosts couple to the gauge field through a derivative coupling 
and this derivative eliminates the constant modes. Now as the ghosts and 
anti-ghosts only occur in the internal loops in the Feynman diagrams, they 
are unphysical fields. So the constant modes will always be eliminated in any 
perturbative calculation. 

We propose therefore to use an effective propagator obtained by first 
subtracting out this constant mode and then taking the zero mass limit of 
the regulated propagator. Note that, we are free to add any finite constant to 
this effective propagator, as the contribution coming from any constant will 
not contribute in the perturbative calculations, due to the reasons mentioned 
above. It appears likely that the use of this effective propagator will lead to 
a consistent perturbative theory. However the consistency of the theory is 
not obvious because it is not clear how the removal of the zero modes from 
the Faddeev-Popov ghosts will affect the BRST symmetry of the theory. 

Thus, we define D^ 2 {x,x') by 

D%M = D m2 (x,x') - — ^. (5.26) 
Then from Eqs. f !5.24p and f !5.26p . we have 

^M=EE y;' )+ t' (5 - 27) 

L=l cr V > 

Clearly this is convergent in the zero mass limit, and its zero mass limit plus 
any arbitrary finite constant C/167T 2 gives us the effective propagator for the 
ghost fields Dq S : 

lim D%(x,x>) + -^~ 2 = Df(x,x>). (5.28) 

This propagator appears for very different physical reasons in the works of 
Allen and Turyn on covariant graviton propagator [6]. 



33 



5.3 Explicit Expression for the Effective Prop- 



agator 

In this section we derive an explicit expression for D$ s using the work of Allen 
and Jacobson [38]. This propagator has been derived and used in a different 
context, as was mentioned above [6]. Following Allen and Jacobson, we define 
fi(x, x) to be the geodesic distance between spacelike separated points x and 
x' in de Sitter spacetime. Also the variable z is defined as, 

z = C os 2 (0 . (5.29) 

In terms of the variable z, the solution to Eq. (I5.19P is given by (see Sec. 2 
ofRef. [SB]) 

D ™< z ) = T^r(a + )r(a_)F[a + , a_; 2; z\. (5.30) 
where a + and a_ are given by 
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a + = 2 + Y4 _m2 ' ( 5 - 31 ) 



a. = (5.32) 

This propagator clearly diverges in the zero mass limit, because T(a_) — > 
oo as m 2 — > 0. However we have seen in the previous section that the 
effective propagator obtained by subtracting the constant L = mode from 
the solution to Eq. (I5.19p . has no IR divergences in the zero mass limit. So to 
find the effective propagator explicitly using the work of Allen and Jacobson 
we shall first verify that the contribution of L = mode subtracted from the 
constant part of the Eq. (15.301) is a finite constant in the zero mass limit. 
Let D c be the constant ^-independent part of D m 2(z) given by Eq. (I5.30p . 

D c = -Lr(a+)r(a_). (5.33) 

l07T z 



We subtract the zero mode contributions given in Eq. (15.231) from Eq. (I5.33P 

D < - 5-TH* = 77h T M T ( a ~) ~ 7T1- ( 5 - 34 ) 
87r z m z lbii z 8n z m z 

We shall verify that this is finite in the zero mass limit. We define 



/; =<\/^-m 2 . (5.35) 
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We have 

We can now write 
r(a + )r(a_) 



P ~2 



2 +P 
\ +P 



-^- + 0(m 4 ), 



(5.36) 



P ) r I 2 ~P 



2~ P 



TT 



P 



P 



Here we have used 



\ +P 



Y{z)Y{l 



2~ P 



shi7r (| — p) 

TT 

— sin 7r (| — p) 
1 



7T 



sin(7T2;) 

Hence r(a + )T(a_) can be approximated as, 

r(a + )r(a_) = ^-^ + (9(m 2 ). 



+ C(m 2 ). (5.37) 



(5.38) 



(5.39) 



The zero mass limit of the constant modes contribution subtracted from D c 
can be written as follows 



lim 

m 2 — s-0 

= lim 



1 



lTT 2 m 2 



m 2 ^0 167T 2 

3 



6 — 3m 2 2 6 
+ C(m 2 ) 



16tt 2 



(5.40) 



Thus this is finite in the zero mass limit and there are no infrared divergences. 
The exact value of this constant is not important as we are free to add any 
constant we like to this effective propagator, because the constant part of the 
propagator does not contribute in perturbative calculations. The important 
point to note here is that it is finite and does not diverge. 

Now we can calculate the effective propagator as follows. We first define 
D' (z) as, 

d 



D'o(z) 



lim 

m 2 ->0 



(5.41) 
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and then define the effective propagator by 



Df(z) 



dzD'(z) + 



C 



16tt 2 



where C is an arbitrary finite constant. Now as 

— F[a, b; c; z] = —F\a + 1, b + 1; c + 1; z] 

dz c 



so Eq. f)5.4ip . becomes 



D' (z) 



1 



16n' 



r r(4)r(l)F[4,l;3;z]. 



(5.42) 



(5.43) 



(5.44) 



With C = 0, Eq. (I5.42p in terms of elementary functions is given by 

1 



16tt 2 



1 . 14 

_-21og(l-,)- T 



(5.45) 



This is the effective propagator that can be used to do perturbative calcula- 
tions. 
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Chapter 6 

Ghosts for Perturbative 
Quantum Gravity in de Sitter 
Spacetime 



In this chapter we shall calculate the ghost propagators for perturbative 
quantum gravity in de Sitter spacetime. We shall again use the fact the 
Feynman propagators in the Euclidean vacuum [I] in de Sitter spacetime can 
be obtained from the corresponding Green's function on the S 4 by analytical 
continuation [55] . 



6.1 Ghosts for Perturbative Quantum Grav- 
ity 



Since the ghosts for perturbative quantum gravity are vector fields, we need 
to review the formalism of Allen and Jacobson for the vector propagators 
(see Sec. 1 of Ref. [38J). Let x and x' be two spacelike separated points and 
let fi(x, x') be the geodesic distance between them as before. One defines 
the unit tangent vectors n a (x, x') at x and n a /(x, x') at x' along the geodesic 
between these two points as follows, 



n a (x,x') 
n a >(x,x') 



V a /i(x,x'), 
V a 'ii(x,x'). 



(6.1) 
(6.2) 



37 



In addition one defines a parallel propagator g ac > such that if A a is a vector 
at x, and if A c the vector at x' obtained by parallelly transporting A a along 
the geodesic, then 

A d = gtA\ (6.3) 
Now as the unit tangents at x and x' point away from each other, so we have 

g a c ,n a = -n d . (6.4) 

One also writes the metric at x and x' as g ac and g a >c', respectively. Now any 
maximally symmetric bi-tensor can be expressed as a linear combination of 
9ac, 9a'c', n a, n c' and g ac i with the coefficient of each term depending only on 
z = cos 2 (/x/2). For example a maximally symmetric bi- vector V ac r(z) may be 
expressed as 

Vac'(z) = a(z)g ad + {3(z)n a n c >. (6.5) 

We have already presented the ghost Lagrangian for perturbative quan- 
tum gravity in general spacetime in Eq. f)3.17p . So we can start from the 
following ghost Lagrangian, 

C 9 h = -iV a c h [V a c b + V b c a -2kg ab V c c c + 
c c V c /i ab + h ac W b c c + h bc S7 a c c - 
kg ab g cd [ceV e h cd + h ce V d c e + h de V c c% (6.6) 

This Lagrangian can now be written in terms of a free ghost Lagrangian 
Cg^ ee and the Lagrangian for interactions The free part of the ghost 

Lagrangian is given by 

4r = -W a c b [V a c b + V b c a - 2kg ab V c c c }. (6.7) 

So the free field equations for the ghosts and anti-ghosts are given by 

V d [V d c a + V a c d -2kg da V c c c ] = 0, (6.8) 
V d [V d c a + V a c d -2kg da V c c c } = 0. (6.9) 

Now we define L a (m 2 ) as follows 

L b a (m 2 ) = V 6 V a - 5 b a a - - m 2 5 b a , (6.10) 

where 

fl = j^T- (6-11) 
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Free field equations, Eqs. ( 16 .8p and (16.91) . in de Sitter spacetime can be 
written as: 



L b a (-6)c b = 0, (6.12) 

L b a (-6)c b = 0, (6.13) 

where 

L>(-6) = V 6 V a - 5 b a D - 2/T 1 V a V 6 + Q5 b a . (6.14) 

Here we have used Eq. f)4.1Up . Here again we choose the Euclidean vacuum 
[I] for calculating the time-ordered product of fields. Let |0) be the Euclidean 
vacuum state, then Feynman propagator for the ghost fields would be given 
by 

(0\T[c b (x)c c ,(x r )]\0) = iG bc ,(x,x'), (6.15) 
where G bc i would satisfy 

L b a (-Q)G bc ,(x,x') = -g ac ,5(x,x'). (6.16) 

Now as in the Euclidean vacuum state the Feynman propagator can be ob- 
tained from Green's function on the four sphere by analytical continuation 
[38] , so to find the Feynman propagator we shall first find Green's function 
on S A . 

On S 4 any smooth vector field can be expressed as a linear combination 
of vector spherical harmonics A^ a and the gradient of the scalar spherical 
harmonics V a Y La (see Sec. 5 of Ref. |5]). The vector spherical harmonics 
satisfy |4"U] , 

-nA? = [L(L + 3) - 
V a A L a u = 0, (6.17) 

where L = l,2,3,4,5-- - and they are normalized as 

J d 4 x^gA La A L,(T ' = S LL '5 aa '. (6.18) 

Here all the quantum numbers on a three-sphere S 3 are denoted by a. The 
degeneracies for these scalar spherical harmonics d s and the vector spherical 
harmonics d v are given by [6] 

d s = i(L + l)(L + 2)(2L + 3), 
6 

d v = ^L(L + 3)(2L + 3). (6.19) 
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The scalar spherical harmonics satisfy Eq. (I5.9P and their gradient covectors 
are normalized as W^ a = V a Y La / \JL{L + 3) by using Eq. (15.91) and Eq. 
(15.101) because 

J d 4 x^g ab W^W b * LW = L{L \ 3) J d A x^gg ab V a Y L °V h Y* L '°' 

= 5 LL '5 aa '. (6.20) 
Now on S* 4 , we have (see Sec. 5 of Ref. [2]) 

g ad 5\x,x') = 6Z)\x,x') + 6i s J\x,x>), (6.21) 

where 

oo 

CAx,x') = EE^M^V), (6-22) 



= EE^w^ iV M- ( 6 - 23 ) 



L=l cr 

OO 



L=l <T 



We can show that this is a valid expression for delta-function by first ex- 
panding any vector function f(x') c > on S 4 in terms of A^r{x') and W^ a {x'), 
then using Eq. (16.181) and Eq. (16.201) . That is, following a similar line of 
argument as was used in the scalar delta-function case, we can show that 

d 4 x'^/g~&jf c '(x')g ac ,6\x,x') = f a (x). (6.24) 

We also have (see Sec. 5 of Ref. [6]) 

G ac ,(x,x') = G v ac ,{x,x') + G s ac> (x,x'), (6.25) 

where 

oo 

L=l a 

oo 

G s a Ax,x>) = ^^W^WH/f'fx'). (6.27) 

L=l a 

Here k\ and k2 are L-dependent constants which can be determined from the 
following equations: 

L a b (-6)Gl,(x,x') = -6 { Z )4 (x } x') } (6.28) 
L a b (-6)G s ac ,(x,x') = -5 { b S /(x,x'). (6.29) 
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Thus we get, 



k 2 



(L + l)(L + 2) -6' 
L(L + 3) 



1 - 



L(L + 3) + 3/3 



(6.30) 
(6.31) 



So we have, 



^ac/{, X i x ') 
Gac'( X i x ) 



EE 

L=l cr 

oo 



A^(x)Af a (x') 
(L + l)(L + 2) -6' 

Y La (x)Y* La (x') 



(6.32) 



L=l cr 

' L(L + 3) + 3/3 



L(L + 3) 



(6.33) 



The right-hand side of Eq. f 1 6 . 3 2 j) is infrared divergent and this infrared 
divergence is caused by the contribution coming from the L — 1 modes 
because (L + 1)(L + 2) - 6 = for L = 1. But if we regulate Eq. (KT2g]l by 
changing m 2 = —6 to m 2 = —6 + p 2 thus adding a small mass p 2 to m 2 , then 
we get 

Lt(-6 + p 2 )G V a y\x,x') = -<£?(*,*')• (6.34) 

By repeating the above procedure we find that the regulated propagator 
G , is given by 



L = l cr 



A^(x)A* c r(x') 
(L + l)(L + 2)-6 + p 2 



(6.35) 



In the zero p 2 limit, G^ p \x, x') — > G^ c ,(x, x') and G^f \x, x') thus diverges 
in the zero p 2 limit. Furthermore, we have seen that this divergence is caused 
by the L = 1 modes. So we can write Eq. (I6.35P as, 



Al°(x)A*}°(x') 



+££ 



^ — — (L + l)(L + 2)-6 + p 2 



2 

Let Q^c' be the solution to (see Sec. 3 of Ref. 

L a b (-6 + p 2 )QC(x,x') = -g M 5\x,x') 



A^( x )Af a (x') 



(6.36) 
(6.37) 
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Then by comparing Eq. (I6.34p with Eq. (16.371) . we can write the solution to 
Eq. flCTj) as (see Sec. 5 of Ref. 0), 



G v a y\x,x') = QC(x,x') + ^^V a V c ,Df(x,x'). (6.38) 

Here we have used Eq. (I6.23P and the fact that L%(— 6 + p 2 ) is invertible. 
Here Dq S (x, x') is the effective propagator given by Eq. (15.281) . as there is no 
contribution coming from the zero modes for D (x,x') due to the action of 
derivatives on it. 

In Eq. (I6.33P there are no zero modes because of the action of derivative 
on the scalar spherical harmonics. So the only IR divergence in Eq. (16.331) 
are caused if we choose (3 = —n(n + 3)/3, where (n = 1,2,3, •••). The 
covariant graviton propagator is IR divergent for the same values of (3 [12J. 
The graviton propagator used by Antoniadis and Mottola [13] corresponds 
to the n = 1 case and is IR divergent as a result. If we avoid these value of 
/3 there will be no IR divergences left in Eq. (16.331) . 

Now because of Eq. (I6.33P we can write the scalar part of the ghost 
propagator as 

G s ac ,(x,x') = -\v a V c \Df{x,x')-D w {x,x')}. (6.39) 
6 

There is no contribution from the L = mode in Eq. (I6.33P because of the 
action of the derivative. So we have used in Eq. (16.391) . the same effective 
propagator that was obtained in Eq. (15.281) for the Yang-Mills case. 



6.2 Effective Propagator 

We have seen in the previous section that the vector part of the ghost propa- 
gator is infrared divergent due to the L = 1 modes. We shall now argue that 
the L = 1 modes do not contribute to the calculations of the time-ordered 
product of the ghost fields. To see that we note that the L = 1 modes are 
the Killing vectors on S A [4"I] . This can be seen as follows. Let a vector f a 
satisfy the Killing equation on S" 4 , 

V a f b + V b f a = 0. (6.40) 

Taking the trace of Eq. (16.401) , we get 

V a f a = 0. (6.41) 
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So f a is a divergencesless vector. Now taking the divergence of Eq. (I6.40p . 
we find 

(-□ - 3) f a = 0. (6.42) 

Which is the equation for L — 1 mode. Conversely, suppose f a satisfies 
Eq. (E22D, then 

J d A x^{V a f b + V b f a ){V a f b + V b f a ) 
= -2 J d A x^f a V b {V a f b + V b f a ) 
= 2 J d 4 x^f a (-a-3)f a 

= 0. (6.43) 

We can conclude from this equation that V a fb + Vb/ a = 0. Thus, (— □ — 
3) f a = implies V a /& + V b f a = 0. In other words, each L = 1 mode is a 
Killing vector and vice versa. 

The interaction part of the ghost Lagrangian can be written as: 

C = -iV a c b [£ c h ab - kg ab g cd £ c h cd ]. (6.44) 

Since the coupling term of the anti-ghosts to the metric perturbation h ab in 
Eq. ( 16.44)) is proportional to V a c 6 + V fe c a , and since the L = 1 modes are 
the Killing vectors, the L = 1 modes do not contribute to the loop diagrams. 
Now as ghosts and anti-ghosts are unphysical fields they only occur in loop 
diagrams and thus there is no contribution from the L = 1 modes in any 
perturbative calculation. 

For this reason, as in the Yang-Mills case, we propose to define the vector 
part of the effective ghost propagator by first subtracting out the contribu- 
tions coming from L = 1 modes from the vector part of the ghost propagator 
and then taking the p 2 — > limit. 

We denote the contribution coming from the L = 1 modes by Q k a p c , (x, x'): 

Q^M^J^ ^T^ - (6-45) 
Now from Eqs. (16. 36p and (16 .45p . we have 

A^ a (x)Af a (x^ 



lim 

p 2 ->0 
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Clearly this is convergent and has no infrared divergence. As the contribution 
V a V c /-DQ ff (x, x') /(6 + p 2 ) had no infrared divergence in the p 2 — > limit, the 
infrared divergences in G^f \x, x') come from Q p ad (x,x'). So we define the 
effective propagator Q e Ei{ Ob . Ob ^ cLS follows: 



Qac'( X i X ) 



lim 

p 2 ->0 



Qac'( X ' X ) Qac'( X ' X ) 



Now the total effective ghost propagator will be given by 



^ac' i X i x ') 



\im\G V a f\x,x') +Gi,(x,x') -Q^(x,x')] 



p 2 -K) 

lim 

p 2 -K) 



Qac' ( X ' X ) Qac' ( X > X ) 

' V a V c ,^ ff (x,a;') 



+ lim 

p 2 ^o 6 + p 1 

1 



-gVaV* [Df(x,x')-D 3 p(x,x')] 
Q^(x,x') + ^V a V cl D^(x,x'). 



(6.47) 



(6.48) 



This propagator can be used to do calculations in perturbative quantum 
gravity. However just like in the Yang-Mills case the consistency of the 
perturbation theory is not obvious as it is not clear how the subtraction of 
the Killing modes will affect the BRST symmetry of the theory. 



6.3 Explicit Expression for the Vector Part 
of the Effective Propagator 

To find the effective ghost propagator in terms of the variable z we can 
proceed as follows: The solution to Eq. (16.371) . in terms of z, is given by (see 
Sec. 3 of Ref. EH) 



Q&(z) = a V{p2 \z)g ad + (5 v ^\z)n a n d , 



where 



a 



v( P 2 ) 



vV) 



-2z(l -z) d 
3 dz 

a y(p 2 )-Y\z). 



2z - 1 



Y (z), 



(6.49) 

(6.50) 
(6.51) 
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Here \z) is given by 



Y\z) = -^^F[b + ^-zl (6.52) 



where 



K = l + J\ + 6-P 2 , (6-53) 



b- = \-\j\ + ^-P 2 - (6-54) 

It may be noted that we have changed the mass from m? = — 6 to m? = 
—6 + p 2 . The propagator given by Eq. (16.491) clearly diverges as p 2 — > 
because then £» — and so T(6_) — > oo. 

The contribution from the Killing vectors in the Allen- Jacobson formalism 
is given by (see Sec. 5 of Ref. [6]) 

QaAz) = jSyi( 2z - l )9aa> + 2{z - l)n a n a ,]. (6.55) 



This can be obtained by substituting 

15 
1Qtt 2 p 



l k = (6-56) 



into Eqs. ( I6.50p and ( I6.5ip . Just like the Yang- Mills case, we first verify that 
7 fc subtracted from the constant part of r y p (z) given in Eq. (I6.52p is a finite 
constant in the zero p 2 limit. Let 7 C be the constant part of r y p (z) given in 
Eq. (I6.52p . then we have 

3r (I + q) T (| - q) 



where q is given by 

We subtract the 7 fc from 7 C : 



^Vgpg, (6 .58) 



3 r(f + g )r(f- g ) is 

7c 7 64vr 2 (6-p 2 ) 167rV 1 ' 

We shall verify that 7 C — 7 fc has a finite zero p 2 limit. Now we have 

q-\ = - P ^ + 0{p% (6.60) 
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We also have 



+ q 



- q 



+ q 



- q 



X 



+ q r 



(6.61) 



We can use the following approximation 



+ q r 



- q 



IT 



7T 



sin7T (i — q) 



sin7r (| — q) 

-^- + 0{p% (6.62) 



Here we have used 



Thus we find 



Hence we obtain 



lim [ 7c - 7 A 

p 2 ->0 



T{z)T(l-z) 



7Y 



sm(Tiz) 



7c = 



64vr 2 p 2 



lim 



60 



87 



p 2 ^0 647T 2 

87 
~320tt 2 ' 



P 2 + ^(P 4 ) 
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(6.63) 
(6.64) 



60 - fp 2 + (9(p 4 ) 



60 
"^2 



(6.65) 



The exact value of this constant is not important because we are free to 
add any finite arbitrary constant to 7 P (z), as the addition of any finite 
constant will only generate Killing contributions which will not contribute 
to the perturbative calculations. The important point to note here is that 
7 P (z) — 7 fc is finite and so there are no infrared divergences. 

Now we can calculate the vector part of the effective ghost propagator by 
first defining 7' (z) as, 



lim 

p 2 ->o dz 



Y\z) 



(6.66) 



and then defining the effective 7 cff by 

1 cG (z)= [ dzi{z) + 
Jo 



C 
64^2' 



(6.67) 
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where C is a arbitrary finite constant. Now as 



4-F[a. b; c; z] = —F\a + 1, b + 1; c + 1: z] 
dz c 



so Eq. (I6.66P now becomes 

Y(*) = 



i r(6)r(i) 



64tt 2 



6 



F[6,l;4;z]. 



(6.68) 



(6.69) 



With C = 11, 7 (z) in terms of elementary functions is given by 



7 eS (^) 



64vr 2 



1 



+ 



L(l-z) 2 (1-z) 



- 121og(l -z) +4 



(6.70) 



The advantage of this choice of C is that it eliminates terms proportional to 
z in a v ^ eS \z) and /3 y ( eff )(z), thus simplifying the calculations. After substi- 
tuting this value of 7 efI into Eqs. f)6.50p and (16.511) . we get 



Qf,0) = a v ^\z)g ad + (5 v ^\z)n a n d , 



(6.71) 



where 



a 



V(eff) 



(*) 



v ^(z) 



64tt 2 
— 12 -f 
1 



64vr 2 
-16 - 



3(1 -z) 2 
4 

3(1-^) 
4 

~3(l-z) 2 
14 



3(1 -z) 



12(2z - 1) log(l - z) 



24(z- l)log(l -z) 



(6.72) 



(6.73) 



6.4 Explicit Expression for the Total Effec- 
tive Propagator 

To calculate the total effective propagator we have to add the contribution 
coming from the scalar part to the effective vector part of the ghost propa- 
gator. The scalar part of the ghost propagator cannot be written in terms 
of elementary functions for a general value of 0. However the correspond- 
ing scalar field equation for /3 = 2/3 is the conformally-coupled massless 
scalar field equation. (An equation of the form [□ — R/6]<f> = is called 
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the conformally-coupled massless scalar field equation as it is invariant un- 
der conformal transformations). For this value of /3 the mass takes the value 
m 2 = 2 and the scalar propagator takes a simple form 

1 1 



Do z) 



16n 2 1 — 2 
Now (see Sec. 1 of Ref. [38]), we have 

-df{z) 



and 



V c ,/0) = -y/z{l-z) 
1 



V„n r / 



2^(1-*) 
So we obtain the following result, 



dz 

{gad + n a n c 



, , 1 df{z) d ( df(z) 

VaVc ' f{z) = 2^z- 9ac ' + il - z) dz[ Z 



dz 



n„n r t . 



Now from Eqs. (16. 74p and (16 .77p we get 
1. 



6 



V a V c ,D 2 (z) = a b (z)g ac , + f3 b (z)n a n c 



where 



a s (z) 



1 



64tt 2 
1 

64tt 2 



3(1 -zf 
2 



+ 



3 1-z 3(1-2 



(6.74) 

(6.75) 
(6.76) 

(6.77) 

(6.78) 

(6.79) 
(6.80) 



So the full effective ghost propagator will be given by the sum of the 
effective vector contributions given by Eq. (16. 78p and the scalar contributions 
given by Eq. (I6.7ip 

Gf c ,{z) = ( a y ^ + a s )g ac , + {P v ^ + P s )n a n c , 
= a eE {z)g ac , + (3 cE {z)n a n dl 



where 



a cS (z) 
(3 eS (z) 



1 



16vr 2 
1 

16tt 2 



1 



3(1-2) 
4 



3(1-2 



- 3 - 3(22 - 1) log(l - z) 
4 + 6(1-2) log(l - 2 



(6.81) 

(6.82) 
(6.83) 



This is the effective ghost propagator for perturbative quantum gravity with 
(3 = 2/3 in de Sitter spacetime, which can be used to do perturbative calcu- 
lations . 
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Chapter 7 

Symplectic Product and 
Wightman Two-Point Function 



In this chapter we shall first discuss canonical quantization of free scalar field 
theory and then generalize these results to higher spin fields. We shall first 
construct the Wightman two-function for free scalar field theory in de Sit- 
ter spacetime by canonical quantization. Then we shall formally see how 
we can use this method to construct the Wightman two-function for higher 
spin fields. From now on two-point function will mean the Wightman two- 
function. 



7.1 Symplectic Product for Scalar Field The- 



We shall review the scalar field theory in de Sitter spacetime in this section 
[39] . However we shall proceed in a way slightly different from [39], so that 
our procedure is easily generalized to higher spin fields. The Lagrangian for 
a minimally-coupled real massive scalar field theory is given by 



ory 



£=^[-V o 0V>-mV], 



(7.1) 



and the classical equation of motion is given by 



(□ - m 2 )(f)(x) = 0. 



(7.2) 
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Now we can define a quantity we call the momentum current ir c as follows: 

, dC 



7T 



-g 



<9V C 



Using Eqs. (17. 2p . we get 



-gm 2 (f)(x), 



(7.3) 



(7.4) 



where 

V c vr c = ^V c [(-g)- l/2 K c }- (7.5) 

If 0, 0' are two solutions of the field equations, and 7r c , ix' c the momentum 
currents conjugate to them, then we define a current J c as follows: 



J c 



-g 



yv c -0V c ]. 



(7.6) 



If dY> a is a surface element of the spacelike hypersurface S, then the sym- 
plectic product on this hypersurface can be defined as follows: 



dZ c J c . 



Using Eq. (I7.4p . we can show that the current J c is conserved, 



(7.7) 



:VJ0V C - 6'Tl* C ] 



iV r \(b*V c (b' -0'V^* 1 



i[(V c 0*.V c 0' - V c 0*.V c 0') + m 2 (0*0' - 0*0')] = 0. (7.8) 



Now if we consider a spacetime region with volume V, bounded by a future 
spacelike hypersurface S + and a past spacelike hypersurface S~, then by 
Gauss theorem we have 



d A Xy/=gV c J c 



and so we have 



dZ r J c 



(7.9) 



(7.10) 



Let us now consider de Sitter metric which was given by Eq. (j4.16p : 

ds 2 = -dt 2 + cosh 2 t[# 2 + sin 2 ^(d6 2 + sin 2 0# 2 )]. (7.11) 
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If we define n c = (1, 0) to be the past pointing unit normal to the constant- 
time spacelike hypersurface £ and 7^ to be the metric on it, then we have 

n c n c = g ab n a n b = -1, (7.12) 

and 

dY, a J a = d 3 x^n a J a 

= c/ 3 xv^j°. (7.13) 
So the symplectic product given in Eq. (17. 7p can now be written as: 

(0, 0') = -i j d 3 x[0*(x, t)7r ,0 (x, t) - 0'(x, t)vr*°(x, *)]. (7.14) 



7.2 Fock Space 



Let {4> n } and {0*} form a complete set of solutions to the field equation 
f)7.2p . and suppose 



(7.15) 
(7.16) 
(7.17) 



The condition given in Eq. (17.161) does not hold in general and so this is a 
requirement on the complete set of solutions to the field equation, Eq. (17. 2p . 
We also chose M nm to have positive eigenvalues only. This again is not always 
true and so this is again a requirement on the complete set of solutions to 
the field equation, Eq. (17.21) . Also in general n,m can be continuous labels. 
However in de Sitter spacetime these are actually discrete labels. 

Now using Eq. (I7.14p the infinite dimensional matrix M nm can be shown 
to be Hermitian: 



% \ rf 3 x[0:(x,t)vr^(x,t)-0 m (x,t)<°(x,t)] 



-i J rf 3 x[0^(x,t)< (x,t)-0 n (x,tX°(x,t)] 



M* 



As is a real field, we can now expand it as follows: 

<f> = ^[a„0„ + a* 



(7.18) 
(7.19) 
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In quantum field theory when <fi and n c are promoted to Hermitian op- 
erators cf) and 7r c , respectively, and the following equal-time commutation 
relations are imposed, 



x,t),7r°(x'^) 
'(x,t),0(x',t) 
[7r°(x,t),7r°(x',t 



z5(x, x') 

0. 
0. 



(7.20) 



then a* and a„ become operators a^ n and a n respectively. Here 5(x, x') is 
defined by 

| rf 3 x5(x,x')/(x) = /(x'), (7.21) 

for any compactly supported smooth function /(x). So we can now express 
as follows: 

} = J2^ + a n ( t ) n]- (7-22) 
n 

For this choice of complete set of solutions to the field equation, Eq. (17. 2p . 
we define a state called the vacuum state |0), as the state that is annihilated 
by a n : 

a n \0) = 0. (7.23) 

It may be noted that here a^ n and a n will not satisfy the standard com- 
mutation relations, however we still call them the creation and annihilation 
operators respectively, in analogy with those for the simple quantum har- 
monic oscillator. Many particle states can be built by repeated action of a^ n 
on the vacuum state. 

It may be noted that the division between {</>„} and {(p^} is n °t unique 
even after imposing conditions given by Eqs. (I7.15p - fl7.17l) [39J. Due to this 
non- uniqueness in division between {4> n } and {</>*}, there is non- uniqueness 
in the definition of the vacuum state also. This can be seen by considering 
{4>' n } and {4>' m } as another complete set of solutions to the field equation, 
Eq. (I7.2p . satisfying conditions given by Eqs. (I7.15I) - (I7.17I) . Now we have 



(7.24) 



Here the vacuum state |0') is the state annihilated by a' n , 

a' n \0')=0. 



(7.25) 
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Now many particle states can be built by repeated action of a'^ on |0'). As 
4> n an d 4>n f° rm a complete set of solutions to the field equation, Eq. (I7.2p . 
we can express <p' n as a linear combination of <p n and 0* , 

= ^[^nm^m + Pnm&m)- ( 7 - 26 ) 
m 

By substituting Eq. 07.261) in Eq. (17.241) and comparing the resulting expres- 
sion with Eq. f)7.22p . we find 

a n = ^ } Qnmfl m ~t~ Pnm a rn\i (7.27) 
m 

a n = ^ ] [ a nm Q m "I" /^nm^m] - (7.28) 
m 

The two Fock spaces based on these choices of complete set of solutions to 
the field equation, Eq. (17. 2p . are different as long as (3 nm ^ 0. In particular 
a n \0') does not vanish because 

in 

of 

but, 

a n \0) = 0. (7.30) 
Thus a n \0') is a one-particle state. In fact we have 

(0'|4a n |0'> = YsYs^nkMrnk- (7.31) 

m k 

We shall use a de Sitter invariant vacuum state called the Euclidean 
vacuum state and its higher-spin analogues in our analysis [4]. The Euclidean 
vacuum is often referred to as the Bunch-Davies vacuum |43j. An advantage 
of using this vacuum state is that the two-point functions reduce to the 
standard Minkowski two-point functions when de Sitter radius is taken to 
infinity, after fixing the geodesic distance between two points [33]. 

7.3 Wight man Two- Point Function for Scalar 
Field Theory 

The Wightman two-point function for scalar field theory is given by 

G(x,x') = (0\4>(x)4>(x')\0). (7.32) 
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This can be written as: 

G(x,x') = J~](Q\(a n <j> n (x) + a^*(z))(q m w (a/) + a] n m (x / )|O) 

n,m 

= ^2<l>n{x)(f>* m {x')C nm , (7.33) 

n,m 

where C nm is the commutator, 

C nm = (0\[a n ,al]\0). (7.34) 

Now as (f) and 7r° are Hermitian, by using the canonical commutation relations 
given in Eq. (I7.20p . we get, 

[(<t>n,4>), (4>,<t>m)] 
= - J rf 3 xrf 3 x'[0:(x,t)7T O (x,t)-0(x,t)<°(x,t), 

0(x',t)7r^(x',t)-0 m (x',t)7r o (x',t)] 
= -ij d 3 xd 3 x'[0;(x, t)7^(x\ t) - m (x', t)<°(x, t)]<J(x, x') 

= -2 J rf 3 x[C(x,t)vr^(x,t)-0 m (x,t)< o (x,t)] 

= M nm . (7.35) 
Now we have 

(<p n , $) = ^ a k M nk- (7.36) 
k 

Using Eqs. (17. 18|) and (17.361) . we also get 

0m) = [{<t>m, 0)] f = £ a\M* ml = 4 M lm- (7-37) 

I I 

Now from Eqs. (I7.35p - (l7.37p . we get 

M nk [a k , a\]M lm = M nm . (7.38) 

kl 

Using Eqs. (1734]) and (17381) . we get 

M nk C kl M lm = M nm . (7.39) 

kl 

This equation in matrix notation is written as, 

MCM = M. (7.40) 
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Assuming that M nm has only positive eigenvalues, i.e., that it is invertible, 

we get 

C = M" 1 . (7.41) 
Therefore the two-point function is given by 

G(x, x') = MxWM)M-l (7.42) 



7.4 Wightman Two-Point Function for Ten- 
sor Fields 

In this section we shall formally generalize what we did for scalar fields to 
general non-interacting higher-spin fields, described by tensor fields. Let 
us denote a tensor field by a shorthand notation, A aia2a3a4 ___ an = Aj. The 
Lagrangian C for this field Aj will be assumed to be a scalar function of 
only Ai and V c Ai, as the inclusion of higher derivative terms will lead to 
problems like non-unitarity of the theory [15]. As we are considering only 
free field theories, we can write the most general Lagrangian for higher-spin 
free field theories as follows: 

C = ~T IcJd V c AiV d Aj - ^S ij AjAj, (7.43) 
where T IcJd and S IJ do not depend on Aj and satisfy 

rplcjd rpjdlc Oj 

S IJ = S JI . (7.45) 
The equation of motion is now given by 

V c [T IcJd V d Aj] - S IJ Aj = 0. (7.46) 
We define a momentum current tt cI conjugate to Aj as follows: 

= ^O^A, = ~^T^ dAj . (7.47) 
From Eq. f)7.46p . we have 

V c vr c/ = -y/=g~S IJ Aj. (7.48) 
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Now if Aj and A\ are two solutions to the field equations, Eq. (I7.46p . then 
we define the current J c as follows: 



J c = --^=[A}ir' cl - A'jir ] . (7.49) 

Now we define a symplectic product on a spacelike hypersurface S c with 
surface element dT, c , as follows: 

(A, A') = J (H1 C J C . (7.50) 

By using Eq. (I7.46p . we can show that J c is conserved, 



V C J C = --^=V c [Ay cl - A'jtt 

= iV c [A*T IcJd V d A'j - A'jT IcJd V ' d A*\ 
= i[V c Aj.T JcJd V d A; - V d A'j.T IcJd V c A* 

+S IJ (A*A'j - A'jA})} 
= 0. (7.51) 

Now if we again consider a spacetime region bounded by a future spacelike 
hypersurface S + and a past spacelike hypersurface E~, then using Gauss 
theorem we can show that 

I dZ c J c = I <fil c J c . (7.52) 

JT.+ JT,- 

Let {A/ n } and {A^ m } be a complete set of solutions to the field equation 
(I7.46p . and suppose 

(A n , A m ) = M nm , (7.53) 
(A n ,A* m ) = 0, (7.54) 
(A* n ,A*J = -M nm . (7.55) 

Here again Eq. (I7.54p is a requirement on the complete set of solutions 
to the field equation f)7.46p . If the theory has gauge symmetry, M nm will 
contain zero eigenvalues and thus will not be invertible [23]. However after a 
suitable gauge-fixing term is added to the classical Lagrangian given in Eq. 
(17.431) . M nm will contain no zero eigenvalues [24J and thus be invertible. In 
this thesis, we will explicitly invert M nm for linearized quantum gravity in 
de Sitter spacetime, in the covariant gauge. 
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It can also be shown that M nm is a Hermitian matrix by repeating the 
argument used in the scalar field theory case. Now after adding a suitable 
gauge-fixing term and imposing canonical commutation relations, we can also 
expand Aj in terms of creation and annihilation operators as follows: 

Aj = Y,KA In + alA* In }. (7.56) 

n 

Here the vacuum state is defined by requiring it to satisfy 

a n |0) = 0. (7.57) 

It may be noted that, just like in the scalar field theory case, there is still 
an ambiguity in the division between {v4/ n } and {v4} n }, and thus there is 
also an ambiguity in the definition of the vacuum state [39] . The Wightman 
two-point function is now given by 

G ir (x,x') = {0\A I (x)Ap(x')\0} 

= J2 A in(x) A i'm(x')(0\K,o>U\0)- (7-58) 

ran 

As after choosing a suitable gauge-fixing term M nm becomes invertible, thus 
following a similar line of argument as was used in the scalar field theory 
case, we can show that 

M-^ = (0\[a n ,al}\0). (7.59) 
Therefore we can now write the two-point function as follows: 

G(x, x') IV = A In {x)A Ilm (x')M-l (7.60) 

nm 

It will also be useful to define a Klein-Gordon type product as follows: 

(A,, A 2 ) = -i J dZ a [A* n V a Al - A I2 V a A{*}. (7.61) 

In the above mentioned scalar field theory case this Klein-Gordon type prod- 
uct is equal to the symplectic product given in Eq. (17.71) . but this is not 
always the case. For example, for a scalar field theory given by the following 
Lagrangian, 

C = -L[-V a 0V> - m 2 2 ], (7.62) 

we have 

(0 1 ,0 2 )=/i(0 1 ,0 2 ). (7.63) 
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Chapter 8 



Spin-2 Field in de Sitter 
Spacetime 



In the previous chapter we derived a formal expression for the Wightman 
two-point function of tensor fields. In this chapter we shall review certain 
properties of linearized quantum gravity in de Sitter spacetime, which will 
be used to explicitly calculate the graviton Wightman two-point in covariant 
gauge in the next chapter. 

8.1 Linearized Quantum Gravity in de Sitter 
Spacetime 

In this section we review the linearized quantum gravity in de Sitter space- 
time. After dropping total divergences, the Lagrangian for linearized quan- 
tum gravity in de Sitter spacetime can be written as [14j, 



This Lagrangian is invariant under the following gauge transformation, 



■gr 




(8.1) 



SAhbc = V c A b + V fe A, 



(8.2) 
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where A c is any vector field on de Sitter spacetime. Now we need to add 
a gauge-fixing term to this Lagrangian for linearized quantum gravity in 
de Sitter spacetime and we choose to add the gauge-fixing term given in Eq. 



So after neglecting total derivatives, the sum of this gauge-fixing term and 
classical Lagrangian for linearized quantum gravity in de Sitter spacetime 
can be written as 



This is the Lagrangian of massless spin-2 field in de Sitter spacetime. It will 
be seen that it is of calculational advantage to first add a mass term and then 
take the zero mass limit of this massive theory at the end of our calculations. 
So during most calculations, we shall start from massive spin-2 field and 
then set the mass to zero at the end of the calculation. The Lagrangian L 
of massive spin-2 field in de Sitter spacetime is obtained by adding a mass 
term called the Fierz-Pauli term, to the massless spin-2 field Lagrangian: 



(USD, 








V a h ac V b h, 



1 P + l 

2 /3a 



be 



) 



xv o w%*-~ h ab h ab + h 2 

+-M 2 [h ab h ah - h 2 ] . 
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The equation of motion obtained from this Lagrangian is 

\v b VaK + v c v a h a b ] 



1 U 11 

-nhhc — 

2 I 2 2a 



03 



a/3" " 2 1 ^ 



(3a 



x[V c V b h + g bc V a V d h 



he + -j9bch 



■-M 2 [h bc - g bc h] 



0. 



5.6) 



8.2 Symplectic Product 

Now to calculate the symplectic product, we first calculate the momentum 
current TT abc , which is given by 



7T 



abc 



dC 



-g 



dV a h bc 



(8.7) 



It can be explicitly written as 

1 



7T 



ahc 



-g 



l V °/i 6c +(1-1.) [g ac V d h db + g ab V d h dc ' 



a^ 2) 9 V V4 2(3a 
: [2g bc V d h ad + g ac V b h + g ab V c h]} . 



Now if h bc and h bc are two solutions to the field equation, Eq. (I8.6p . and if 
TT abc and ir' abc are the momentum current conjugate to them, then we find 
using Eq. (I7.49p that the conserved current J a is given by 



J a 



-fj 



:[h cb 7T - tl cb 7T \. 



19) 



The symplectic product defined on a spacelike hypersurface S is given by 



(MO = / dz a r. 



(8.10) 



Now because of Eq. (17.511) . we have 



V r J c = 0. 



(8.11) 
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So for a spacetime region bounded by a future spacelike hypersurface S + and 
a past spacelike hypersurface £~, because of Eq. ( I7.52p . we have 

d£ c J c = / dZ c J c . (8.12) 



8.3 Mode Decomposition 



In this section we shall decompose h bc into modes. In the next chapter 
we shall calculate the graviton Wightman two-point in covariant gauge by 
calculating the contribution to it coming from each of these modes. 
We can decompose the field h bc into scalar, vector and tensor parts 

hbc = Cbc + Abe + E bc . 
Here the scalar part C bc is given by, 

Cbc 



(V b V c - -g bc u)B + ^g bc h. 



(8.13) 
(8.14) 



The vector part, Ab c is given by 

A bc = V c A b + V b A c , 

and satisfies, 

A b b = 2V b A b = 0. 

The tensor part E bc satisfies 



El 



V C E, 



be 



0. 

= 0. 



(8.15) 
(8.16) 

(8.17) 
(8.18) 

Now we want to find the equations of motion for each of these parts. We 
start from the scalar part C ab . First by substituting h ab = C a b, where C a b is 
defined in Eq. (18.141) . into the equation of motion (18.61) and then using the 



(8.19) 

(8.20) 

(8.21) 

(8.22) 
(8.23) 



formula for R bac given in Eq. (I4.10p . we find the following equation, 



g^X.h + X 2 B] + VaVbfrB + Y 2 h] = 0, 



where 



x 2 



(/3 + 1) 2 + 1 



a/3 2 



4a/3 



3 

°-4 



3M 2 



1 _ 3(1 + 0) 
4 

3 — a 



4a/3 



□ 2 + 3 
M 2 



a(3 



M 2 

□ + -^-n, 



4a 
a/3 -3/3 

4a7T 



□ + 

a 

4 
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Now if we define 



X 1 h + X 2 B = 
Y\B + Y 2 h = 

then we can write Eq. (18.191) as 

Qabli + V a V b z/ 
Taking the trace of Eq. (18.261) . we obtain 

1 



0. 



-Uv. 



Now its divergence, we obtain 

V«[n + 4]z/ = 0. 

So we have 

[□ + 4]zy = 4k, 
where k is a constant. Thus we get 

v = fi + k. 

Now a particular solution to Eqs. (18.241) and (|S2 

h = h {p) = 0, 

2a 



can be written as 



B 



B (p) 



-k + 



2/i 



8.24) 
8.25) 

8.26) 
8.27) 

8.28) 
8.29) 

8.30) 

8.31) 
8.32) 



6-M 2 a 2-M 2 ' 

Any solution Eqs. ( I8.24p and f )8.25j) can be written as a sum of this particular 
solution and the general solution. The general solution is the solution to the 
homogeneous equations obtained from Eqs. ( I8.24p and (I8.25P by setting 
/i = v — 0. However, C a b = for the particular solution h = h 1 ^ and 
Since we take the M 2 — > limit in the end, the special cases 
M 2 = 2 and M 2 = 6/a are not relevant here. Hence we only need to 
consider the solution to the homogeneous equations given by 



X 1 h + X 2 B 
Y\B + Y 2 h 



Thus we have 

1 _ 3(1 + 0) 
4 



+ 



4a/3 

(/3 + 1) 2 



.4 a/3 



af5 2 
3 — a 



Aa 



3 

□ + - 

a 



P + l 

4a(3 
M 2 



l 

" 4 
B + 



□ 



0. 
0. 



M 2 ' 
□ + — □ 

3 3M 2 
a/3-3/3-4 



(8.33) 
(8.34) 



B 



Aa(3 



0, (8.35) 
0. (8.36) 
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After eliminating the fourth-order derivative of B, we can rearrange these 
equations as: 

UB = a 1 B + a 2 h, (8.37) 

ah = a 3 h + a 4 B, (8.38) 



where 



«2 



12 - 2aM 2 

a — 3 

a/3-3/3-4 

£(a-3) ' 



6(a/3 - 3/3 - 2) / M 2 
° 3 a -3 V ~ ~2~ 

-9/3(a/3-3/3-4) / M 2 \ / « 2 \ 

The equations of motion for the vector part are obtained by substituting 
h a b = A a b, where A a b is defined by Eqs. ( I8.15P and f)8.16p . into the equations 
of motion f )8.6p . Thus we get, 

V 6 K + V o y 6 = 0, (8.40) 

where, 

V a = V b [V b A a + V a A b ] - aM 2 A a . (8.41) 
The general solution to this equation is 

V b [V b A a + V a A b ] - aM 2 A a = f a (8.42) 

where f a is a Killing vector and thus satisfies 

VJ6 + V b f a = 0. (8.43) 

So a particular solution is given by 

Af = -Jfc (8.44) 

However A ab = for this solution, A a = Aa ■ So following a similar line of 
argument as was used for the scalar part, we can set f a = 0. Thus we get 

V b [V b A a + V a A b ] - aM 2 A a = 0. (8.45) 
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Now by using Eq. ()4.1Up . we can write this equation as 



[□ + 3 - aM 2 ]A a = 0. (8.46) 

The equation for the tensor part is obtained by substituting h ab = E ab , where 
E ab is defined by Eqs. (I8.17P and (18.181) . into the equation of motion (I8.6p . 
Thus we get 

[□ - 2 - M 2 }E ab = 0. (8.47) 



8.4 Harmonic Modes on de Sitter Spacetime 

To find the contribution of the vector and tensor parts to the two-point 
function we expand the vector and tensor parts into harmonic modes on 
de Sitter spacetime. The spatial part of these modes is given by the spherical 
harmonics on S 3 and its temporal part is given in terms of the associated 
Legendre function. We denote the coordinates on S 3 by x and define \ using 
Eq. P~T71) as x = tt/2 - it. 

Now we review some properties of scalar, vector and tensor spherical 
harmonics on S 3 [4"0] . Let the indices on de Sitter spacetime be denoted 
by (a, b, c - • • ) and on S 3 be denoted by (i, j, k ■ ■ ■ ). We denote the scalar, 
vector and tensor spherical harmonics on S 3 by Y 1(7 (x) , Yf a (x) and Yf? (x) , 
respectively. Here a denotes all the indices other than £. Thus a denotes 
all the quantum numbers on a two-sphere S 2 from now on. It may be noted 
that a previously denoted all the quantum numbers on a three-sphere (in 
the chapters on ghost propagators). We also denote the covariant derivative 
on 5* 3 by Vj, the metric on S 3 by rjij and VjV* by V 2 . For scalar spherical 
harmonics on S 3 , we have i = 0, 1, 2, 3 ■ • • and 

- V 2 F te = £(£ + 2)Y ea . (8.48) 

For vector spherical harmonics on S* 3 , we have I = 1, 2, 3, 4 • ■ ■ and 

-V 2 Yf a = [e(£ + 2)-l}Yf°, 

VV/ CT = 0. (8.49) 

For tensor spherical harmonics on S 3 , we have £ — 2,3,4,5 ■■ ■ and 

[£(£ + 2)-2}Yff, 
0, 

0. (8.50) 



V 2 Y l ° 

V ij 
7 1 ij 
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Next we review some properties of the associated Legendre functions P v 
The degree v can be lowered or raised as 



1 -x d 



d 



i d 

dx 



+ vx ) P 



{l-x')—-{u+l)x)P l 



(x) = {v- /j,)P v ^(x), 

(x) = _(^ + M+ i)p- + M( x ) 



Now we define D„ as follows: 



d 



D n = — + n cot x- 

dx 



Then we have 

- sinxA^(cosx) 



;i - cos 2 x)- 



d 



— n cos x 



dcosx 

So now because of Eqs. (I8.5ip and (18.53 p . we have 



-sinxD u+1 p-» 



[v - n)P v » x {x), 

-( 1 y + l , + l)p-^(x). 



We also have 



because 



D m sin" xPy M (cos x) = sin" xD m + n P u M (cos X ), 



d 



dx 
sin n x 



+ m cot x 
d 



sin" xi^cosx) 



dx 



+ (m + n) cot x 



P-^(cosx). 



(8.51) 



(8.52) 



P-^(cosx). (8.53) 



(8.54) 
(8.55) 



(8.56) 



There are two types of solutions to Eq. (18.461) denoted by A£ £cr with n = 
0, 1. These modes have £ = 1, 2, 3 • ■ ■ . There are three types of solutions to 
Eq. (g37D denoted by Eg* with n = 0, 1, 2. These modes have i = 2, 3, 4 • • • 
(see Sec. 2, 4, 5, 6 of Ref [H]). Here a denotes all the quantum numbers on 
a two-sphere. We let the solutions to the minimally-coupled massive scalar 
field equation in de Sitter spacetime be denoted by S ta with £ = 0, 1, 2, 3, • • • , 
where 

[□ - m 2 ]S la = 0. (8.57) 

We shall write the solutions on de Sitter spacetime in analogy with the 
solutions on S A , exploiting the similarity of de Sitter metric given in Eq. 
(I4.16P to the metric on S* 4 , after defining x by Eq. (I4.17P as x — 71 / 2 ~ tt- 
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So first we review the defining properties of scalar, vector and tensor 
spherical harmonics on S A (50] ■ For scalar spherical harmonics on S 4 , we 
have 

- uS Ua = [L(L + 3)]S L£a , (8.58) 

with L = 0, 1, 2, • ■ ■ and £ < L. For vector spherical harmonics on S 4 , we 
have 

-UA^ = [L{L + 3) - 1]A^, 

with L = 1,2,3, •• • and £ < L. For tensor spherical harmonics on S A , we 
have 

-DEg* = [1(1 + 3) - 2}El e ^ 
V b E^ a = 0, 

g bc E™° = 0, (8.60) 

with L = 1, 2, 3 ■ ■ • and £<L. 

What we use is the fact that the positive-frequency mode function for the 
Euclidean vacuum in de Sitter spacetime are obtained by solving the same 
equations as these spherical harmonics under the transformation x — Tc/2—it, 
if we allow L, L and L to be non-integers and unconstrained by the condition 
L,L,L > £ (see Ref. [17]). Comparing Eq. QgagJ with Eq. flgTSg} , Eq. 
(JEHD with Ec l- ASH, Eq. f lSTBTj) with Eq. jjg3g]l and using the similarity 
of de Sitter spacetime with 5* 4 , we can write L, L and L as, 

1(1 + 3) = — M 2 , 
L(L + 3) = —aM 2 + 4, 

Z(L + 3) = -m 2 . (8.61) 
Now L, L, L to the first order in M 2 or m 2 , are given by 

M 2 

L « - — , (8.62) 
aM 2 

L w 1-^—, (8.63) 
5 

2 

L w (8.64) 

In the limit M 2 — > 0, we have L = and L = 1, for de Sitter spacetime. 
Now the scalar, vector and tensor harmonic modes in de Sitter spacetime are 
given by (see Sec. 2, 4, 5, 6 of Ref. [47]). 
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4 fa ( X ,x; M 
^( X ,x;ilP 

£?x(x,x;M 
£ 2 f( X ,x;M : 



Pf CT ( X ,x;M : 
4^( X ,x;M 2 
4f( X ,x;M 2 



P^f ( X ,x;ilP 

P^( X ,x;M 2 
P°f( X ,x;M< 
P<f (x,x; M 2 



7! 



(sin x) P 



ip-(W), 



cos X )F to (x), 



L+1 



0. 



P; 



-(*+!), 



L+1 



cosx)^(x), 



sin X )- 2 P L X J (cosx)^(x), 



DiP-f+^lcos^V^^x), 



+ 2) 



(8.65) 
(8.66) 

(8.67) 



0. 
0. 



(8.68) 



(smx)- 1 P T l e i +l \cosx)Y l ^( 
1 



(i-l)(i + 3) SinxD2P ^ +1){C0SX) 
x(V,lf (x)+V^(x)), 

L+l 



(8.69) 



(sin X )- 3 P T ^ +1) (cos X )^(x), 



sin ^-^iPf^ (cos X ) V^x) 



£(£ + 2) v ^ 1 L +! 

(^ + fe)^(x), 
3 



f--(sinx) 1 



2(£ + 3)(£-l) 

+ J(iT2) sinxD2Dl ) p i+i +1) ( cos *)' 

_! ' +1) (cosx), 



-p: 



3sinx L+1 
^ ^ + 2) 



(8.70) 



We emphasize again that we have merely used the analogy between 
de Sitter spacetime and S" 4 and not really analytically continued from S 4 
to de Sitter spacetime. If for the modes given by Eqs. fl8.65l) - fl8.70l) . we let 



0,1,2,3,4-.. ,L 



1,2,3,4,5-.. ,L 



2,3,4,5,6- •• and L,L,L > £, 



where I is a positive integer, then they would be regular at x — 71 an d thus 
Eqs. fl8.65p - fl8.70p would define valid spherical harmonics. However, in case 
of de Sitter spacetime these modes do not satisfy all of these conditions, 
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and as a result they are singular at \ = n - But, this point corresponds to 
imaginary time t = iir/2, and hence it is not physically relevant in de Sitter 
spacetime and so we can safely use Eqs. fl8.65l) - fl8.70l) for calculating graviton 
Wightman two-point function in de Sitter spacetime. 

These modes given by Eqs. fl8.65l) - fl8.70l) . regarded as functions on de Sit- 
ter spacetime, can be used to decompose the scalar C ab , vector A ab and tensor 
E ab parts of the field h ab in terms of mode functions. The scalar, vector and 
tensor mode functions that constitute the fields C ab) A ab and E ab , respec- 
tively, are orthogonal to one another with respect to the symplectic product 
as long as M ^ [48J. This implies that C ab) A ab and E ab commute with 
one another upon quantization. Thus the quantization of each of these fields 
can be considered separately. However in the limit M — > 0, the vector modes 
constituting A ab , and tensor modes constituting E ab , with n — 0, 1 coincide 
[TO] . For this reason, we shall consider the quantization of the vector and 
tensor parts, A ab and E ab together, but consider the quantization of the scalar 
part C ab , separately from the vector and tensor parts. 
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Chapter 9 



Wightman Graviton Two-Point 
Function 



In this chapter we shall explicitly calculate the contributions to the covariant 
Wightman graviton two-point function coming from the scalar and vector- 
tensor sectors. 



9.1 Physical Equivalence 

In linearized quantum gravity the Wightman two-point function of the gravi- 
ton field h a b has no physical meaning by itself because this theory has gauge 
invariance under the gauge transformation, 

S A h bc = V c A b + V b A c , (9.1) 

where A c is any vector field on de Sitter spacetime. One can find tensor fields 
at x that are linear in h a b and are invariant under this gauge transformations 
[IE] . An example of such a tensor field is the linearized Weyl tensor W a b c d(x) 
given by 

W abcd (x) = W [ab][cd] (x), (9.2) 
W ab cd{x) = V c V b h a d(x) + h cb (x)g ad (x). (9.3) 

Now, if a Wightman graviton two-point function Gb c b'c'{x, %') can be writ- 
ten as 

G bc b'c'{x,x') = P bc b'c'{x,x') + Qbcb'c'{x,x'), (9.4) 
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where 



Qaba'b' = VaKL'b* ( X > X> ) + V &-^aa'&' ( X > X ') 

+Va>K 2 M (x,x') + V b ,K 2 aba ,(x,x'), (9.5) 

for some Kl a , b ,(x,x') and K 2 bb ,(x, x'), then the two-point function of a local 
gauge-invariant tensor field linear in h bc will be the same whether one uses 
Gbcb'c'(x,x') or P bcb ' c /(x,x') as the graviton two-point function. 

This motivates the following definition: We say that the two graviton 
two-point functions, G bcb ' c '(x,x r ) and P bcb ' c '(x,x') are physically equivalent 
in linearized gravity if they differ from each other by a two-point function of 
pure gauge form Q bcb > c i(x } x'). 



9.2 Symplectic Product for Vector and Ten- 
sor Modes 

The Klein-Gordon type product given in Eq. f)7.6ip . for vector and tensor 
harmonic modes on de Sitter spacetime is given by (see Sec. 8, 9 of Ref. [II]), 

4L(Z + 3)<5"> CT ' 



(E U °,E U °) 
(E oe °,E oi '°) 
(E 2i °,E 2ew ) 
(A oe °,A oiW ) 



(£ -!)(£ + 3)T(£ - L)T(£ + 1 + 3)' 
3 [1(1 + 3) + 2] 1(1 + 3)5"> ff ' 
'£ + 2)(£- 1)(£ + 3)T(£ - L)T(£ + 1 + 3)' 

5 ll> 5 a*> 



T{£- L)T{£ + L + 3) 

2(1(1 + 3) + 2)£"> g ' 
~ £(£ + 2)T(£ - L)T(£ + 3 + L)' 



The Klein-Gordon type product (E ne ° ,E n ' e ' a ') and (A nttT ,A n ' e ' <T ') vanish for 
n 7^ n'. The vector mode functions corresponding to A™ ta are given by 

A^ = V a Af° + V b A n a ia . (9.7) 

As shown in Appendix A, the symplectic product for the vector and tensor 
modes is related to their Klein-Gordon type product as follows: 

(A nea ,A n ' e '' J ') = -M 2 {A nta ,A n ' l ' (T '), (9.8) 
(£*,£ n '' v ) = ^(E n£(T ,E n ' e ' a '). (9.9) 
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Now recall 



L 
L 



1 - 



aM 2 



M 2 



(9.10) 
(9.11) 



It will turn out that we need L and L only to the first order in M 2 . Using 
Eqs. flEBD and (KTE]) . we get 



(A oe °,A oi '°') 
(A U °,A U '°') 



2M 2 5"> CT ' 



(£ -1)(£ + 3)T(£ + M 2 /3)T(£ + 3 - M 2 /3) : 

r2//\ / f2 n\ sll' £<t<7' 



3M 2 (M 2 - 2)5" 5° 



2£(£ + 2)(£ - 1)(£ + 3)T(£ + M 2 /3)T(£ + 3- M 2 /3) ' 

rll' roV 



5" 5° 



T(£ + M 2 /3)r(£ + 3 - M 2 /3) ' 

2M 2 (aM 2 -6)S ll 'S arT ' 
~£(£ + 2)T(£ - 1 + aM 2 /5)T(£ + 4 - «M 2 /5) ' 
2M 2 <5"> fJ ' 



(9.12) 



r(£ - 1 + aM 2 /5)T(£ + 4 - aM 2 /5) ' 
To simplify these expression further, we note 

T(£ - M 2 /3)T(£ + 3 + M 2 /3) 
= (£ + 2 - M 2 /3) (£ + 1 - M 2 /3) (£ - M 2 /3) 

xr(£-M 2 /3)r(£ + M 2 /3), (9.13) 

and since T(£ — M 2 /3)T(£ + M 2 /3) is analytic and even in M 2 , we have 



r(i - M 2 /3)r(£ + m 2 /3) « [r(£)] 2 + c(m 4 ) 

Hence, to the first order in M 2 , we find 

T(£ + M 2 /3)r(£ + 3-M 2 /3) 



(9.14) 



(£ + 2)\{£-l)\ 



1 



^ 1 \M 2 
^ P + k) 3 



,fe=0 



(9.15) 



Similarly, we have 



T(£ - 1 + «M 2 /5)r(f + 4 - «M 2 /5) 
+ 3)!(£-2)! 1- ( 



aM 2 



(9.16) 
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We have found the symplectic product for vector-tensor harmonic modes 
in de Sitter spacetime, in this section. We have also seen that the Wightman 
two-point function for any free tensor field can be obtained from this sym- 
plectic product by using Eq. (17.601) . So in the next section we shall use the 
symplectic product calculated here to find the vector-tensor contribution to 
the graviton Wightman two-point function. 



9.3 Vector- Tensor Contributions 

Now we can write the vector-tensor contribution to the Wightman two-point 
function by using Eq. (I7.60p as follows: 

^ r aba'b'\ X i- L ) tj^niu £nlo\ 

n=0 i a ^ ' ' 

I f VV g^)gVU 
2.^1 2-^ 2-^ (firda fine.o\ ' \ J - 1 ') 

n=0 t cr V ' ' 

where for (n = 0, 1), we have defined A^ a as, 

A n a ^ = S7 a Af° + V h Af°. (9.18) 

In Ref. [TU| a physical Wightman graviton two-point function P a ba'b>{%i x') 
was obtained. To do so first de Donder gauge V b h ab — V a h/2 = was 
imposed. Then it was shown that the trace could be gauged away and thus 
the conditions h a a = and V a h ab = were also imposed. Finally it was 
shown that all the modes apart from E^{x) could be gauged away. Thus 
it was shown that only E^ b (x) modes contribute to the physical part of the 
Wightman graviton two-point function. So P a ba'b'(x, x') thus obtained was 
given by 



We will show here that this physical Wightman graviton two-point function 
given in Ref. [10] is physically equivalent to the Wightman graviton two- 
point function in covariant gauge. To do so, we write down the tensor- vector 
part of the Wightman graviton two-point function in covariant gauge, 

G[h] a = V b h ab - (1 + /T 1 ) VA (9.20) 



as 



G abZl'( X i X> ) = P aba'b'{x,x') + M aba , V {x,x') + L aba > b > (x , x') , (9.21) 
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where, 

1 4< 



. 1 . . _ Anla ( \ A*nl<i ( l\ 



^afea'fe' (S, X') = 



n=0 cr 

1 x Ante I ~\ A*nAa(~f\ 



n=0 £=2 cr ^ ' 2 

, ^afc l X W / q9 oN 

n=0 £=2 cr V ' ' 



Now substituting Eq. (19TT2]) into Eq. (19T23]) . we get 

, v^v- 1)(^ + 3)„ ( n M 2 \/ M 2 \ 



=2 cr 

1 „ /„ aM 2 \ „ / aM 2 

r + — — r ^ + 4 — — 



2M 2 V 5 

2£(£ + 2)(£ -!)(£ + 3) 



3(2- M 2 )M 2 



+EE 

^=2 cr 

xr + x) r (' + 3 " x) ^to^SfV) 

£(£ + 2) / «M 2 \ /„ aM : 

+ ^tt^ K - 1 + — — r U + 4 - 



2M 2 (6-aM 2 ) V 5 / V 5 

xAT(x)A*f/(x')] , (9.24) 

The contribution from each mode to L a &a'6' is divergent in the zero mass 
limit. However we shall see that the terms of order 1/M 2 cancel out. Now 
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for small M 2 , substituting Eq. (KTT5]1 and Eq. flgJ5) into Eq. 1$FZty . we get 

i °° r / 1 1 2 1 \ 

£=2 o- L V fc=o / 



1 00 

+T9EE^ + 2 )( £ + 3 ) ! ^-2) ! 



=2 (T 



V fc=o / 



. (9.25) 



It may be noted that if E^ u = /2(£ - 1) and E 1 ^ = A 1 ^" /(£ - 1), in 
the zero mass limit then the 1/M 2 singularity will be absent in this limit. In 
fact we show in Appendix B that to the first order in M 2 , we have 



Elite 



2{£-l) 

1 

(F=T) 



(9.26) 



The exact form of fc™^ which is M 2 -independent, for n — 0, 1 is also given 
in Appendix B. Now if we substitute Eq. (19. 26j) into Eq. (I9.25p . and define 
Ar a t%,forn = 0,las 

Tvrnta I l\ _ \nl<j I \ A*nlo l l\ , Ante/ u*tifa/ A , ynla I \ a *nte ( l\ 
aba'b'K i ) ab \ X ) A a'V \~ J ^ ab \ x )^a'b' \ X ) ' K ab W^Vb' \ X > ' 

(9.27) 

then in the M 2 — >■ limit, we have 

00 

Wvfox') = 53J3 [w oe N°J:, b ,(x,x') +w u N£, v (x,x')] , (9.28) 

1=2 IT 
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where 



w e 



12 



£{£ + 2)(£ + Z)\(£ - 2)\ x 



1 1 x ^ 1 a a x ^ 1 

~9 + q TXT + fi" _ r" 2^ 7T 



3 + fc 

fc=0 



1 



(£ + 3)!(£-2)! 



5 ^ £ + k 

fc=— i 



^ £ + fc 3 ^ £ + k 

k=-l k=0 



(9.29) 



We can now write L aba i h i as, 

Laba'b' (X, x') = VaVla'V ( X , X> ) + V 6^aa'6' ^0 



+V 't&y(z, z') + V 6 '^ 6a ,(a:, O, 



(9.30) 



where for n = 1,2, 



1 oo 



n=0 £=2 a 
1 oo 



4#oo + 



(9.31) 



n=0 ^=2 ct 

There is a finite contribution to the Wightman graviton two-point func- 
tion coming from the £ — 1 vector modes. In the small mass approximation 
this finite contribution can be written as 



M aba , v (x,x') 



X 



a^(x)a:^(x') 

2M 2 



X 



o /tow \ /i*oict/ 
2M 2 (6-aM 2 ) 



aM 4 



(9.32) 



nlcr 



n=0 cr 

-60. This has a finite zero mass limit because 
V a A b llcr + Vf,A" lCT = in the zero mass limit for £ — 1, since for 

nlcr 



where 6 = 30 and fei 



L = 1 these are Killing vectors. Expanding A™^ , we get 

dA r ^ u (x) 



Af°{x)=Af°{x) \ L=l + {L-l) 



dL 



\L=l ■ 



(9.33) 
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Now from Eqs. fIDTTUj) and (JjT53)» . we get 

= < CT (x) | L=1 + (9.34) 

where 

Thus ^4^ CT i s °f or der aM 2 . So in the zero mass limit, we obtain 

1 

M aba , v {x,x') = ^^ab n (V a X? 1 ° + V^ 1 °)(V al X2 1 °* + V v Xp t >*) 

— V a x^ a , b / (x, x ) + Vf,rr ao / 6 , (x, x ) 

+V a 'i(^0 + V&'xL'O^ A (9.36) 

where 

sum = EE^^ nlfT ( v «' x ^* + v ^"' 1<T *)' 

n=0 a 
n=0 <t 

Thus, we have calculated the vector-tensor contribution to the Wightman 
graviton two-point function in this section. In particular, we have shown that 
it is physically equivalent in linearized gravity to the physical Wightman two- 
point function of Ref. [10J. In the next section we shall calculate the scalar 
contribution to the Wightman graviton two-point function and find that it 
is pure gauge. 



9.4 Equations of Motion for Scalar Modes 

The scalar contribution to the Wightman two-point function has been cal- 
culated for a(3 — 3/3 — 4 = in Ref. [TJ]. The scalar contribution to the 
Wightman graviton two-point function was found to be a pure gauge contri- 
bution for these value of a and 0. We shall calculate the scalar contribution 
to the graviton Wightman two-point function for arbitrary values of a and 
/?. 

Let us recall that the field equations for scalar modes, given by Eqs. (18.371) 
and (I8.38p . are 

UB = a 1 B + a 2 h, (9.38) 
Uh = a 3 h + a 4 B, (9.39) 
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where 

12-2aM 2 , , , 

ai = 5— , (9-40) 

a/3-3/3-4 

° 2 = fla-3) ' (9 - 41) 

6M-3/3-2) / M 2 \ 
% = ^ I 1 "";' (9 ' 42) 

a 4 = zMM^ij U -*£)(!- » . (9.43) 
a — 3 \ 2 J \ 6 / 

In deriving these equations we have assumed that a ^ 3, and /5 ^ 0. 
However, the final result for the scalar contribution to the Wightman graviton 
two-point function is well defined even for these values of a and /3. 

Eqs. ( I9.38P and ( I9.39P suggest that B and h are non-trivial linear com- 
binations of two scalar fields with different masses unless a/3 — 3/3 — 4 = 0, 
which leads to = 0. Let a linear combination of B and h, h + XB, 

satisfy the free scalar field equation with a definite mass. Now we have 

u(h + XB) = a 3 h + a 4 B + X(a x B + a 2 h) 

= (a 3 + Xa 2 )h + (a 4 + Xa x )B. (9.44) 

As we want Eq. f !9.44j) to satisfy a free scalar field equation, with a definite 
mass, the right hand side of Eq. (I9.44p must be proportional to h + XB. That 

is, 

^±P± = (9-45) 
(a 3 + Aa 2 ) 

Thus we get the following quadratic equation in A, 

A 2 - ^^A - ^ = 0. (9.46) 
a 2 a 2 

Let Ai and A2 be the two roots of this equation. Then the product of the 
roots is given by 

AiA 2 = (9.47) 
a 2 

and the sum of the roots is given by 

A a + A 2 = (9.48) 
a 2 

Now the scalar fields <j>i and (fi 2 defined by 

01 = h + XtB, (9.49) 

2 = h + X 2 B, (9.50) 
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satisfy 



□0i = (a 3 + Aia 2 )0i, (9.51) 
□02 = (a 3 + A 2 a 2 )0 2 - (9-52) 



Solving Eqs. (I9.49[) and (I9.50P for B and h, we get 



B = f*-^, (9.53) 

A 2 — Al 

* = A ^~^ . (9.54) 
A 2 — Ai 



9.5 Scalar Contributions 

Since the scalar part C ab {x) depends on two scalar fields B{x) and h(x), we 
need to find (0\h(x)h{x')\0) , (0\B(x)B(x')\0) and (0\B(x)h{x')\0) and take 
the limit M ->■ in order to find (0|(7 a6 (x)(7 a / 6 /(x / )|0). 

Let two solutions to the equations of motion Eq. (18.61) be Cf£ and C^, , 
given by 

Cl: = (V b V c -\g bc n)B^ + ^g bc h^, 
C(% = {VvV c >-\gv c iU')B l '°' + ~^g Vd h e °'. (9.55) 

Here the covariant derivatives V a and V / only act on quantities at x and 
x', respectively. We also define □' = V a /V a ' and □ = V a V a . Also let 
7Tg C ^ abc (x) and Kg u ' a b c {%') be momentum current conjugate to and 
Gy^, , respectively. We can now calculate the scalar contribution to the sym- 
plectic product by substituting these in Eq. (18. 9p . This symplectic product 
for the scalar part is calculated in Appendix A, and is given by 

[C l % C £V ) = -K x {B la , B l '°') - K 2 {h £a , h e ' a '}, (9.56) 

where 

K ' = At 3 " M2 )( 6 " aM2 )' ( 9 ' 57 ) 

4(a — 6) 

and where (B ea ,B e,r7 ') and (h ea ,h ta ') are defined by Eq. fl73T|) . We note 
that B ea and h £r7 are linear combinations of modes with different masses, so 
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the Klein-Gordon-type product (B ea , B ia ) and (h ea , h l ' a ') are not conserved 
separately. 

Now let and be positive frequency modes satisfying Eqs. (I9.5ip 
and (I9.52p . respectively. Let us also require that 



\V1 7 01 CT ) 

then by using Eq. (I7.42p . we have 

A m? (*,x') = ^^0f(x)^(x'), 



(9.59) 



A m 2(s,a/) = J^X/ 



>f(x)<^(x') 



(9.60) 



Here A m 2(x,x') is the standard two-point function of the scalar field with 
squared mass m 2 , (see Sec. 2 of Ref. [3B]) and 



777,1 = a 3 + A x a 2 , 



"^2 = a 3 + ^20-2- 

Now we can expand <f>i(x) and <j>-i(x) as 



x a 



02 (x) = ^^0f(x)a 2te + 0^(x). 



£ a 

We now write the field C a f,(x) as follows: 

C ab (x) = C 1 ab (x) + C 2 ab (x), 

where 



with the differential operators x\ b and £ 2 6 defined by 



-1 



x 



ab 



X 



ab 



X 2 — Ai 
1 

A2 — Ai 



V a V fe - 7^6 □) - TfafcAl 



(9.61) 



(9.62) 



(9.63) 



(9.64) 



(9.65) 
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We now expand the fields C\ b {x) and C 2 b ( 



C 2 ab (x) 



X a 



*2for/ \ f 



x a 



2&r> 



(9.66) 



£ cr 



where 



C a b a ( X ) — X ab4>l{ x )i 

Now from Eq. (19. 56p . we get 

K x + \lK 2 
(A 2 - Ai) 2 ^ 
K x + A?AT 2 



(C u °,C 2e '°') 



(A 2 - Ai) 



0. 



(9.67) 

(9.68) 

(9.69) 
(9.70) 



Now we note from Eqs. (EH]), (jg3SD and (19T4TL 



— AiA2-^2 = —-^2 
«2 

3 



4(a - 3) 



(2- M 2 )(6-aM 2 



(9.71) 



Substituting Eq. fl9T7Tj) into Eqs. f l938|) -f l938|) .and using Eq. fl939|) . we get 

(C 1& ,C W ) 



A2-^2 ,r<7'\ 
1 )H / 



(C 2fo ,C 2 « v ) 



Ai — A2 

A2-^2 fiM'fiaa' 
Ai — A2 
A l-fG 1 ,1a ,ot„'. 



A2 — Ai 



A 1-^2 ^«y ff ' 



A2 — Ai 

Thus the only non-zero commutators are given by 

Ai — A2 



(9.72) 



X2K2 

A2 — Ai 
AiK 2 



OfpO„„i. 



(9.73) 



SO 



Now using Eqs. (17.601) and ()9.73p . we get 

(0\C ab (x)C a , b ,(x')\0) = 



l cr 



X2K2 



■ ("ilia 1 \^i*Uu 1 t\ 



+^J 1 C%°(x)C^(x>) 



(9.74) 



Now from Eqs. (ICTj) . (1931) and (l9~74j) . we get 

(0\C ab (x)C a , v {x')\0) 

= EE 



l cr 



a 2j ft 2 

Thus we have 



A 2 X 2 

'^ab^Vb' ^mj (x, X ) + 



Al^ 2 



' X ab4 > 2 a { X ) X a'b'^2 ° ( X 



X ' " " * X lb X a'b'^m?;( X i x ) 



(oi^x^xOlo) 

(O|0 2 (x)0 2 (x')|O> 



AiiT 2 

Ai — A 2 

A 2 /v, 
A 2 — Ai 

Ai^ 2 



(9.75) 



A m? (x,x'), 



A m2 (x,x'). 



(9.76) 



Next we use these to find the Wightman functions for B(x) and h(x). We 
define A av (x, x') and A dlff (x,x') as: 



A av (x,x') 
A dlff (x,x') 





(x 


x') 


+ A m 2 


(x, 


x') 








2 








(x 


x') 




(x, 


x') 






to to 




— mf 






A m 2 


(x 


X 7 ) 


- A m ? 


(x, 


x') 



(9.77) 



a 2 (A 2 - Ai) 

Then using Eqs. (19. 53 p . ( I9.54p and ( 19. 76 p . we obtain 



(9.78) 



(0\h(x)h(x')\0) 

(0\B(x)h(x')\0) 
(0\B(x)B(x')\0) 



^[A av (x,x') 

J^2 

V ~A dm (x,x') 



—a 2 



Ai + A 2 



K 2 
1 

K2 



fl2 A diff (x,x'), 



(9.79) 
(9.80) 



AiA 



-A av (x,x') 



(9.81) 
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We have, from Eqs. (ICTjl . ( jg37D and (ESHJ) 

1 



= -/3>-3), (9.82) 
^ = -p( a p-3p-4). (9.83) 



In the zero mass limit, we obtain from Eqs f l9.40p - fl9.43p 



M 2 -X) 2 

This implies 



lim AiA 2 = 9/T, (9.84) 

Af 2 ->0 

l im Al + Aa = -3/3, (9.85) 



lim Ai = lim A 2 = -30, (9.86) 

M 2 ->0 M 2 -s.O 



and hence by Eq. (19.611) . we get 



lim m? = lim ml = 313. (9.87) 

A/ 2 ^0 Af 2 ^0 



Now we define A l(x, x') as 



Note that 



Af 2 ^0 

So in the zero mass limit, we have 



A^(x,x') = -^A m ,(x,x'). (9.1 



lim A av (x,x') = A 3/3 (x,x'), 

M 2 ^0 

timA^&x') = A^{x,x'). (9.89) 



(0\h(x)h(x')\0) = f3 2 (a-3)A w (x,x') 

-3f3 2 (af3 -3/3 -4)A ( ^(x,x'), (9.90) 
(Q\B{x)h{x')\<d) = -/3(a/3-3/3-4)Aj ) (x,x'), (9.91) 

(0\B(x)B(x')\0) = -^-A 3 p(x,x') 



- aP -f- A A^x>). (9.92) 

By differentiating the formula (□ — m 2 )A m 2(x,x') = with respect to m 2 , 
we obtain 

nA^ 2 (x,x') =m 2 A^ 2 (x,x')-A m 2(x,x'). (9.93) 
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^aba'b'i X i X ) ~~ ^aba'b'( X i X ) + ^aba'b'( X i X ) + ^abJb'( X i X ) 

where G^^x 7 ), G^^x') and G^) v (x, x') are given by 



(9.94) 



So the scalar part of the graviton Wightman two-point function, which is 

G { 2'A X , X 1 = (0\C ab (x)CM x 'M, 

can be written as 

<S) / n _ MSO) , n , MSI) 



(9.95) 



G (S0) , n _ 
(j aba>b'\ x i x I — 



GibJb'( X ' X ) 



Va'Vf/ - -g a 'b'{x)W 



1 



x [V a V b --g ab (x)nj (0\B(x)B(x')\0) 
lV a V b V a ,V b/ (0\B(x)B(x')\0) 



+ ^{9ab^a'^b' + ga'b'^aSb) 

+PW -3/3-4) Ag(x, *') 



-A 3/3 (x,x') 



-Y^9ab9a'b'[-f3{af3 - 3/3 - 8)A 3/3 (x,x') 
+3/3 2 (a/3-3/3-4)Ag(x ! :r')], 



^ (^V a V 6 - -<7 ab nJ (0|£(aO/i(aO|0) 
-I/3(a/3 - 3/3 - 4)^ 6 ,V a V 6 Ag(x,x') 



+— p (a /3 - 3/3 - 4)s ab5a , 6 , [-A 3/3 (x, x') 
+3/3Ag(x,x')" 



G%2vM = T ^9ab(x)9a'b'(x f )(0\h(x)h(x')\0) 

= YQ9ab(x)g a 'b'(x') [/3 2 (a - 3)A 3/9 (x,x') 
-3/3 2 (a/3-3/3-4)Ag(x,x') 



(9.96) 



(9.97) 



(9.98) 



(9.99) 



The contributions proportional to g ab g a > b > cancel out and the scalar part of 
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the Wightman graviton two-point function can now be written as 
d&vM = V a V 6 V a ,V 6 ,(0|5(a;)S(x')|0) 

= ^ 'a S la'b' ( x , x ) + VbS l aa i b i (x, x') 

+V a ,s 2 M (x,x , ) + Vvs 2 M (x,x'), (9.100) 

where 

sl, b ,(x,x') = ±V b V a ,V v (0\B(x)B{x')\0) 

+75 , a'fe'VfeA 3 ^(x,a;'), 
o 

s 2 aW (x,x') = ~V a V b V a ,(0\B(x)B(x')\Q) 

+lg ab V b ,A 3 p(x,x'). (9.101) 
o 

Thus, we have calculated the scalar contribution to the Wightman gravi- 
ton two-point function and found it to be of pure gauge form. 



9.6 Full Graviton Two-Point Function 

Now the full Wightman graviton two-point function can be written as: 

Gaba'b' {x,x) = G a b a 'b' ( X i X ) + ^ikt'b' i X l X ) 
Paba' 

b'{X,X) + V a S ba 'b' {X, X ) + VbS aa , h t [X, X ) 
+Va'4^,l') + V 6 '4, a ,(x,x') + V a vl a , v {x,x') 
+Vbvl„,, V (x, x') + V a iV 2 aW (x, x') + V b 'vlba'( x , x> ) 
+^aXba'b' { X , X ') + ^bX aa i b ' (x, X ) + V a 'X abb , ( X , x ') 
+Vb' x lba> (X,X') 
= Paba'b' (X, X') + Qaba'b' (x, x') , (9. 102) 

where Q a ba'b'(x, x') is a pure gauge contribution given by Eq. (19. 5p with 

K ba'b> = S ba'b'+ V ba'b'+ X ba'b' aIld K lbb> = S lbb> + V lbb> + X lbb' ■ The physical tWO- 

point function obtained by fixing all the gauge degrees of freedom in Ref. [TU] 
is denoted by Pb c b> c >{ x , x ')- Therefore the covariant Wightman graviton two- 
point function is physically equivalent in linearized gravity to the physical 
Wightman two-point function. 
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We established the equivalence of the vector-tensor part of the Wightman 
graviton two-point function to the physical two-point function P aba i b i(x,x') 
of Ref. [10] in detail using the mode expansion. In Ref. [17] it was shown 
that the scalar sector is of pure gauge form for a specific value of a and (3. 
Here we have generalized this to arbitrary value of a and (3. 

Now suppose that a covariant Wightman graviton two-point function has 
an infrared cut-off e and that it is divergent in the limit e — > 0. As it is 
physically equivalent in linearized gravity to P a ba'b'(x, x'), and as P a ba'b'(x, x') 
is not infrared divergent, the two-point function of a local gauge- invariant 
tensor field will not depend on e, i.e. will not be infrared divergent. This 
also means that the infrared divergences in the Wightman graviton two-point 
function for a certain gauge choice, for example as noted in Ref. [13], will 
not appear in the Wightman two-point function of any local gauge-invariant 
tensor field. So in this sense these infrared divergences in the Wightman 
graviton two-point function are gauge artefacts in the context of linearized 
gravity. 
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Chapter 10 
Conclusion 



This thesis consisted of two parts. In the first part we studied the Feynman 
ghost propagators for Faddeev-Popov ghosts for Yang- Mills theories and per- 
turbative quantum gravity in de Sitter spacetime and in the second part we 
analyzed the equivalence of the covariant Wightman graviton two-point func- 
tion to the physical one obtained in Ref. [TUj . 

We found that the ghosts and anti-ghosts in Yang-Mills theories satisfy 
the minimally-coupled massless scalar field equation and as a result the Feyn- 
man propagator for Yang-Mills theories in de Sitter spacetime is IR diver- 
gent. We also found that the Feynman propagator for perturbative quantum 
gravity is also IR divergent. 

However it was shown that if we regularize these propagators by adding 
a small mass then the modes responsible for these IR divergences do not 
contribute to the time-ordered product of the fields. We thus found the 
effective propagators by first subtracting these regularized modes which cause 
the IR divergences and then taking the zero mass limit. This way we obtained 
the effective IR finite propagators for Yang-Mills theories and perturbative 
quantum gravity. 

It seems likely that the use of these effective propagators will lead to 
consistent perturbative theories. However the consistency of these theories is 
not obvious because it is not clear how the removal of some of the modes from 
the Faddeev-Popov ghosts will affect the BRST symmetry of these theories. 
We can not naively remove the constant or Killing modes from these theories 
because if we do so, they will be brought back as BRST transformation of 
some other modes. This point needs further investigation and we would 
have to modify the original BRST transformations for consistency. It will 
be interesting to find the modified BRST transformations for these theories 



after constant or the Killing modes are removed. 

We have dealt only with integer-spins and it will also be interesting to 
derive the propagators for the half-integer spin ghost fields in supergravity. 
At present it is not clear if these propagators will be IR finite or IR divergent. 
It is also not clear that if they are IR divergent, then we can still use the 
method we used for integer spin fields to remove the IR divergences. 

We analyzed the Wightman graviton two-point function with arbitrary 
gauge parameters in the second part of this thesis. We found that the co- 
variant Wightman graviton two-point function is physically equivalent to the 
physical one obtained in Ref. [10], in the sense that they result in the same 
two-point function of any local gauge-invariant quantity. We then argued 
that as the physical Wightman graviton two-point function obtained in Ref. 
[TU] is known to be IR finite and the covariant Wightman graviton two-point 
function is physically equivalent to it, so any IR divergences in the Wightman 
graviton two-point function, for example as noted in Ref. [13], has to be a 
gauge artefacts at least in the context of linearized gravity. 

We should emphasize that we have only studied the free theory and at 
present have nothing to say about Wightman two-point functions in the 
interacting theories. It will be interesting to analyze what happens when 
interactions are introduced. 

It will also be interesting to analyze the equivalence of the covariant 
Wightman graviton two-point function to the physical one in anti-de Sitter 
spacetime in arbitrary dimensions. In fact we can also analyze this equiv- 
alence for Wightman graviton two-point functions in four dimensional flat 
spacetime. This is widely believed to hold but has not been done yet and 
it will be interesting to see if a result similar to the one obtained in case of 
de Sitter spacetime holds for flat spacetime also. 



87 



Appendix A 
Symplectic Product 



A.l General Procedure 

We want to relate the symplectic product for scalar and vector modes to the 
Klein-Gordon type product for them. For this purpose the following formula 
is useful 

V a V fe 0Vy - (0 <-> 0') = -V a 0D0' + (0 <-> <j>'), (A.l) 

where (0 -H- 0') is defined to mean that the preceding expression is repeated 
but with <f> replaced by 0' and 0' replaced by 0, and = indicates equality up 
to terms that vanish when they are integrated over the space. To show this 
we first note that 

V a V 6 0V b 0' - (0 <-> 00 = Q a ~ [V a 0D0' - (0 <-»• 0')], (A.2) 

where Q a is a total divergence of an antisymmetric tensor: 

Q a = V b [V a 0VV - V 6 0V a 0']- (A.3) 

This total divergence of an antisymmetric tensor can be converted into a 
surface term by Gauss's divergence theorem on the spatial section S 3 because 

J d^ b V a F ab = J rf 3 xv^V a F a0 

= j d^di{^gF i0 ). (A.4) 

As there is no boundary on S 3 , we just drop such terms. We have used Eq. 
( jZH5D here. 
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A. 2 Symplectic Product for the Scalar Modes 



We shall calculate the contribution to the symplectic product coming from 
scalar sector here. Let B £cr ,B ecr ,h £cr ,h ecr be solutions to Eqs. (18.37'j) and 
flOBD and let 



n i'a' 
^bc 



1 

4' 



■^Qbch 
- A 9b c ri + -B bc 



(A.5) 



where 



B, 



B 



fa 
be 
it 



be 



{AV b V c -g bc n)B"% 
(4V b V c - g hc u) £? £V 



(A.6) 

We denote Cfc, B[° c , B iu , h l ° by C ab , B bc , B, h and C 6 f Bjf <J ', 5* £V , fr* €V 
by C' ab , B' bc B', h' here for simplicity of notation. However it may be noted 
that they are modes and not fields which will be promoted to operators. 

Momentum current conjugate to C bc can be found by substituting h cb = 
C bc in Eq. (18. 8p . Now if the momentum current conjugate to C bc is tt?^, 
then we have 

(A.7) 



7T 



abc 
(C) 



_.ofec i „abc 
(B) + *7 



where we define to be the contribution coming from the pure trace part 
of the scalar modes which is obtained by substituting h bc = Qbch/A in Eq. 
(18. 8p and as the contribution coming from the traceless part of the scalar 
modes which is obtained by substituting h bc = B bc /4 in Eq. (18.81) . Thus we 

get 

4 + 3/3 - a/3 



„„abc 



„—abc 
~9 (B) 



8a/3 



1 1 + /3 
r 1 4 4a/3 



(g ab V c h + g ac V b h) 
(1 + /3) 2 



2 



a/3 2 



V C V°B - -g cb uB 
4 



1 1 



(A.f 



X 



V a V°5 - -o d6 nS 
4 



V a V c 5 - -q^uB 
4 



+ 



1 10 + 1 
"2 /3a 



V d V a fi - -g da uB 
4 



We can simplify by using 



V d V a 



-g da a)B = -V a (a + A)B, 



(A.9) 



(A.10) 
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as follows, 



r ,„abc 

-g^{B) 



1-^ ) (g ab V c + g ac V , )(a + A)B 



1 1 + 8 

+- I -o + — f 

2 ap 



3 
4 

~V { V n V - -g*u j D 



g bc V a (a + 4)5 



(A.ll) 



If Jq is the contribution to J a defined in Eq. (18. 9p . from the scalar modes, 
then we have 



where 



ra 

J c 



labc s~if „abc 

;[O c6 7r (c) - G cb 7T 



ja 

J (l) 



ja 

J (2) 



ja 



1 



1 



[hg C b^[ a B) - h'gcb^(B) 



labc t}I „abc 



[Bbc^B) ~ B , 



6c 7r (B)J 



(A.12) 
(A.13) 

(A. 14) 
(A.15) 
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Thus we have 



4 + 3/3 - a(3 1 1+j (1 +/j) 2 ' 
16o/3 + 4 + 1^/3 a/3 2 
x{hV a ti -h'V a h), (A.16) 



J ( a 2) = 4 + ^ , fl ^ (y a V b BV b ti - W a W b B'W b h) 



4a/3 

± + W-^ {nBV a h ,_ nB , v a h) 



16a(3 

34 + 3 ^- a VvC + 4) B 



16 a/3 

-W a (D + 4)5'), (A.17) 
Jf 3) = l^ 1 -^) (V a V 6 5- V 6 (D + 4) J B' 



-V a V h B' ■ V b (D + 4)fi) 

-A fl- -1 (□BV a (D + 4)5 / 
16 \ ay 

-□S'V°(n + 4)B) 
-i(V 6 V c 5 • V a V b V c B' - V b V c B' ■ V a V b V c B) 

+i(nfi • V a nB' - uB' ■ V a uB). 

(A.18) 

To simplify the term W b W c BW a W b W c B' - (B B') in J ( a 3) , we first note that 

V b V c B ■ V a V b V c B' 
= V b V c B ■ [(g ca V b B' - g cb V a B') + V b V a V c B'] 
= V a V b B -V b B' -nBV a B' 

+V b V c B ■ V b V a V c B'. (A.19) 

Then by considering this term together with the term obtained by inter- 
changing B and B', we have 

V b V c BV b V a V c B' - V b V a V c BV b V c B' 
= V b {V b V c BV a V c B' - V a V c BV b V c B') 

-nV b BV a V b B' + V a V b BnV b B' 
= -V b uBV a V b B' + V a V b BV b uB' 

-3V b BV a V b B' + 3V°V b £V b B'. (A.20) 
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From Eqs. (IA.19|) and (|A,20p . we have 



V 6 V c £V a V 6 V c £' - (B -H- 5') 
V a V b BV b B' - nBV a B' - V b nBV a V b B' 
+3V a V b BV b B' -(B<* B') 
-nBV a B' + V a V b BV b {n + A)B' - (B £') 



So we can now write as 



J ( a 3) _(!__] V°V 6 SV 6 (D + 4)5' 



"^f 1 "-! aBV a (n + 4)B' 
lb \ a J 

+ loBV a (n + A)B' B') 
3 



a 



1 

16 



J V a V & 5V fe (n + 4)5' 
J DBV(D+4)B' 



-(5 <R> 5'). 
Now we note that using Eq. (I8.37p . we have 

□ + 4 5 = - 2 - M )B + ^ % 

a — 3 p(a — 3) 



(A.21) 



(A.22) 



(A.23) 



So the terms of the form V a V b B'V b B' or V a Vb-BV b /i' or these with the 
primed and unprimed functions swapped in Jfy and are 

4 + 3/3 - a/3 (vavfe5V ^, _ vav&5 , v ^ ) 



4a/3 



+t ( i - - ] V fc V a 5V 6 



3 
a 

2-M 2 



V b V a B'V b 



J^(2-^+ a/3 - 3/3 "V 
a-3 V ; /3(a-3) 



(V a V h BV h B' - V a V b B'V b B). 



a — 3 



/3(a-3) 



(A.24) 



Now from Eq. flA.ll) . we have 



V a V b BV h B' - W a W b B'W a B = -V a BuB' + V a B'nB, 



(A.25) 
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so we can write 



where 



ja 

J (hh) 



ja 

J (BB) 



ja 

J (Bh) 



ja I ja rsj jQi | j a , ja 
J (2) "+" J (3) — J (BB) "+" J (Bh) "+" J (M)' 

3 (a/3-3/3-4) 2 



16 a/3 2 (a - 3) 
3(2-M 2 )(6-aM 2 ) 
4 a - 3 

4 + 3/3 - a/3 12 - 2aM 2 



16a/3 a — 3 

3 4 + 3/3 - a/3 2a 



{hV a ti -h'V a h), 
(BV a B' - B'V a B), 
(BV a ti - B'V a h) 



(A.26) 

(A.27) 
(A.28) 



16 a/3 a — 3 
2 - M 2 a/3 - 3/3 - 4 



{2 - M 2 ){h'V a B - hV a B') 



/3(a-3) 



(-V a Bti + V a 5'/i) 



1 

16 

2a 



1 - 



a 



12 - 2aM 2 a/3-3/3-4 
a — 3 /3(a — 3) 



(BV a ti - BV a h) 



Now we have 



ja _ 
■J(Bh) ~ 



—(2 - M 2 ) Q/3 3/3 4 (W^ - fc'V'B) 
a — 3 p(a — 3j 



4 + 3/3 - a/3 6 - aM 2 la-3a-3/3-4 



(A.29) 



8a/3 a — 3 
x(BV a ti - B'V a h) 

2 - M 2 a/3 - 3/3 - 4 3 4 + 3/3 - a/3 



+ 



+ 



2 /3(a-3) 8' /3(a-3) 
1 a - 3 a(2 - M 2 )(a/3 - 3/3 - 4) 



a /3(a — 3) 2 

(2 - M 2 ) 



(6 - aM 2 ) 



8 a (a - 3)/3 

x(hV a B' -h'V a B) 
= 0. 



(A.30) 



So the terms cancel out. The total B-B contribution is given by J^ BB ^ 
and the total h-h contribution is given by J? hh \ + J?u, which is given by 



J (i) + J (hh) 



1 

16 



(3/3 + 4) 2 3(a/3 - 3/3 - Af 



a/3 2 



a/3 2 (a -3) 



1 



>(a-3) 



(hV a ti -h'V a h). 



(hV a ti - /A7 a /i) 
(A.31) 
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So if 

C ab = + fv a V 6 - ^g ab n J B , 
C' ab = \ti +(v a V b -\g ab njB' , 

then we can write, 

(C,C) = -K 1 (B,B')-K 2 (h,h'), (A.32) 

where 

K x = ir ^-r(2 - M 2 )(6 - «M 2 ), (A.33) 
4(a — 3) 

* = -^T3)' < A ' 34 > 
This is the symplectic product between two scalar modes. 



A. 3 Symplectic Product for the Vector and 
Tensor Modes 

The equation of motion for the vector modes, Eq. (18.46p . can be written as 

[□ + 3 - aM 2 }A b = 0. (A.35) 

Let A^ la and A b 1 ' e ' a ' be two solutions to this equation. We will denote A^ ta 
by A b and A™ 1 a by A' b for simplicity of notation. However it may be noted 
that they are modes and not fields which will be promoted to operators. 
The momentum current conjugate to A bc = V c A b + V b A c , can be obtained 
by substituting h bc = A ac in Eq. (18.81) . Now if the momentum current 
conjugate to A bc is 7T?Jn, then we have 



—abc —abc i —abc IK qc\ 

(A) = ^(A4) (A5) ' {A.3b) 
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where, 



„^abc 



■-V a (y b A c + V c A b ) 



d A b + W b A d ) 

0&V7 fX7'd AC I V7 C A<i\ 



+-[g ac V d (V d A b + V b A 



+g ab V d (y d A c + W c A d )\ 
\[g bc V d (y a A d + W d A a )}, 



„„abc 



— [g ac V d {V d A b + V b A d ) 
2a 

+g ab W d iy d A c + W c A d )] 
' l \g bc V d (y a A d + V d A a )}. 



(3a 



(A.37) 



(A.38) 
(A.39) 



If J a A is the contribution to J a defined in Eq. ( I8.9p . from the vector modes, 
then we have 



ja 



^=[A cb Tr'ffi — A' cb 7T(A)] 



where 



ja 

J (4) 



ja 

J (5) 



--[A bc 7l 



labc 
(A4) 



A I „abc 1 
/i bc 7r (A4)J ) 



-g 



--[A bc TT 



labc a I „abc 



(A5) 



l bc n (A5)\- 



(A.40) 

(A.41) 
(A.42) 



It is well known that the symplectic product vanishes for h ab = A ab in the 
theory without the gauge- fixing term (see Ref. [M]). This implies that 



(A.43) 



- % I dZ a J ( a 4) = 



So the only contribution to the symplectic product comes from J%y Now we 
can simplify y/—g^ A2 ) by using Eqs. (IA.35P and f )4.10p as 



„„abc 



-—[g ac (D + 3)A b + g a \u + 3)A C 
2a 



(3a 



[g bc (D + 3)A a ] 



M 



(3 + 1 



_ {g ac A b + g ab A c ) + >^J_ g bc A a 



(A.44) 
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Substituting Eq. (IA.44p into Eq. (jA.42j) . we get 

J ( a 5) = -M 2 [A b V a A' b - A' b V a A b ] 
Thus we have from Eqs. (IA.43P and f)A.45|) . 



(A.45) 



'(5) 



iM 2 J dL a [A b V a A' b - A' b V a A b ] . ( A.46) 



So we have 



{A, A') = -M 2 {A,A'). (A.47) 
We can write the equation of motion for tensor modes, Eq. (I8.47p . as 



□ 



M 2 ]E, 



ab 



0. 



(A.48) 



nlcr 



Let E b l a and E b ^ a be two solutions to this equation. We will denote E bc 
by E bc and E^ ° by E' bc for simplicity of notation. However it may be noted 
again that they are modes and not fields which will be promoted to operators. 
The momentum current conjugate to E bc , can be obtained by substituting 
h bc = E ac in Eq. (18.81) . Now if the momentum current conjugate to E bc is 
7r"J?j, then we have 



„^abc 



-V a E bc . 



(A.49) 



So we have 



J a E = --[E bc V a E' bc - E' bc V a E 



1 



6cl 



If Jg is the contribution to J a defined in Eq. (18. 9p . from the tensor modes, 
then we have 

(A.50) 
(A.51) 



:e,e') = -(e,e'). 



A. 4 Summary 



In summary, 



—M 2 (A nia , A n ' e ' a '), 
— (^E n ^ cr E n '^' a '^ 



(A.52) 
(A.53) 
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and if 

c£ = \g a ^ + (v a v 6 - \g ab ^j B?° , 

(ft = \ 9ab h^ + {v a V b -\g ab U^B^\ 

then we can write, 

{C l % C e ' a ') = -K X {B^, B e '°) - K 2 (h e °, h ew ), (A.54) 

where 

Ki = 3 „, (2-M 2 )(6-«M 2 ), (A.55) 
4 (a — 6) 

K 2 = l — 

2 /3 2 (a-3) 
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Appendix B 



Explicit Form for the Tensor 



B.l General Procedure 

We shall calculate explicitly k^ a , for n — 0, 1 here, for completeness. It may 
however, be noted that the explicitly form of k r ^ a is not necessary for the 
main results of this thesis. 
Recalling that 

A n a l° = V a Af° + V b A n a e °, (B.l) 

we have 



fcX , for n = 0, 1 




= 2D A^, 
A nl ° = D_ 2 Af° + ViA 



^ = ViAf CT + VjAf ° 



(B.2) 
(B.3) 
(B.4) 
(B.5) 
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So we have 



A 



Ha 

XX 



AW 



A 



Ua 

Ua 
ij 



D.oP, 



2-fL+l 



COS X )^(x 



P L ii +1) (cosx) 



A 



Oicr 
XX 



3 \~> • '3*1 
-2 n d-(W)/ \v<» 



2(sin X )- l D_ 2 p-^>{cos x)F to (x) 



(B.6) 
(B.7) 
(B.8) 
(B.9) 



O^cr 



A Ma 



+ 2) 



1 



sin x\ 



P^ipatxW'W, (B.10) 



+ 2 



+2 cot X P4i +1) (cos x)%'^ ct (x). 



(B.11) 



From Eqs. fl8.66l) - fl8.70l) and (1B.6j) - (1B.llj) we see that the S 3 tensor structure 
for E%£ c is similar to A^ a for n = 0, 1. 

''aft 



Now is defined as: 



K ab 



M- 2 [2{l-l)Ef£-Ag] 
M- 2 [(£ - 1)E% ~ ^] ■ 



We first recall that from Eq. (I8.52p . we have 

d 

D n = — + n cot x, 
dx 



(B.12) 
(B.13) 



(B.14) 



and from Eq. (18. 54p . we have 

-sinx^-.P^ 
-sin xDv+iP' 1 * 



-(i/ + Ai + l)P^(x). 



(B.15) 



B.2 Explicit Form for the Tensor k\[ a 



We can now calculate k}?° as follows: 

AA 



C^^^-iK?-^] = ' (B.16) 
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To calculate klf 7 , we note that 



Xi 

-(M-i), 



2- r L+l 

sinx 



cos x) 



D_ {l+1) P-1[ +1 \cos X ) + (L— l)cot X p-|r^(cos X ) 



sin x 
[(L-l)-(€-l 



-(<+!), 



cos x) + (£ - 1) cot xP L +i +1) (cos x) 



smx 



i 



smx 
L - 1 



-P_.\ y (cosx) H — : 



L+l 



smx 



P L (m) (cosx)-P x ^ +1) (cosx) 



+ 



smx L 
t-\ 



cos xP L +i l) ( cos X) - P L ( cos X) 



(M-l), 



smx 



L+l 



qM 2 



cosx) - ~ r~ {(€- 1)^M 2 |-P-^(cosx) 
smx [ 15 az/ 



P x (m) (cosx) -cosxP 2 ( ^ +1) (cosx) 



From this we find 



Ha 



where 



Qi(x) 



I 



£-1 
sinx 



u=l 

(B.17) 
(B.18) 



(^-l)^^ (m) (cosx) 



MT" ' 15 &/ 



!/=l 



Pr { ' +1) (cos x) - cos xiY (m) (cos x)] } • (B. 19) 



Similarly to find k}j a , we note that 



-(<+!), 

L+l 

sinx 



P» 2 pr^7^(cosx) 



-P 



L+2)^W) - ^cotxP x X lj ( cos x) 



(M-i)/ 



smx 



3) 



-(M-i) 



smx 



L+2 



(cos x)~ L cot X^r+i +1) ( cos X) 



L+l 



smx 



P-^(cosx) + - 



smx 



-(<+!), 



L+2 



cosx) 



-(<+!), 



L+l 



cos x) 



cos X^x+i +1) ( cos X) ~ P t+o'' ( cos X) 



smx 



smx 

P^T'WsxH 



L+2 1 



3 „_(^+i) , 1 



(1 + 3)%^M> |-P-(^)(cosx) 
smx I 15 <7z/ 



M 2 



P 2 (m) (cosx) - cosxPi ( ^ +1) (cosx) 



^=2 

(B.20) 
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Hence we find 



1 



ViY} l ° + V.F 



(lforO)' 



3 i 



(B.21) 



where 



Q2 (x) 



3« — 5 d 



15 <9z/ 
1 

~3(£ + 3) 



p-( m )(cos X ) 



i/=2 



P 2 - ( ' +1) (cos X ) - cosxPr i£+ij (cos X ) • (B.22) 



B.3 Explicit Form for the Tensor k®i a 



Next to find , we use Eq. (|B.17p . and get 



"xx 



Then we have 



1 

M 2 

— ^ 1) 
sin % I 15 <7z/ 



[2(^-1)^-4?] 
3a — 5 9 



— P-^(cosx) 



i/=i 



2(<-l) 



Pi 



•(<+!), 



cosx) — cos xP 2 ^ +1 ' ) (cos 



■X)]} 



F 



. * gi(x)^ 

sin x 



Of a 



(B.23) 



where 



„0&r 
"Xi 



M 2 



+ 2) 



lim — — - 

AP->0 M 2 



X(x)V,F 



1 



D_ 2 D X + 



+ 2) 



sin 2 x 



P, 



,-(<+!) 
-(<+!) 



2(^-1)— AP^r^(coBx) 
sm x i+i 



L+l 



Noting that 



qi{x)= lim 



1 sin x 
APi- 1 



1 „_ 



smx 



P^ +1) (cosx) - P- 2 P L 7r iJ (cosx) 



cos x) 



-(<+!)/ 



(B.24) 



(B.25) 



(B.26) 
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we find 



2 -^-^M>D iqi {x) 
smx 

+ cotx — + 6- v 
dx dx 



M 2 X( X ) 
r d 2 



sin 2 x 



^ii +1) (cos X ) 



d 2 d lT iWr n , (£+l) 2 

+ [6-(L + l)(L + 2)]P L ^ +1) (cos X ) 



P L --J! +1) (cos X ) 



aM 2 P 2 - (m) (cos X ), 



in the M — > limit. Hence 

^iD igi ( X ) + aP-^(cos X ). 



Finally to find /c° £o ", we write it as 

z(x) 



where 



Z(X) 



W( X ) 



k 0£a _ 1 

ij ~£ + 3 



lim — — ■ 

M 2 ^0 M 2 



£(£ + 2) 



Y 



0£a 



(B.27) 
(B.28) 

(B.29) 



sinx 



3P 2 Pi - 



+ 2) 



sin 2 x 



xp _(, +1)(cosx) _ 2{£ + 3)jDi p L -(m) (cosx) 
P 2 P 1 -^ + 2 )~ 2 



(B.30) 



lim — — ■ 

M 2 ^0 M 2 



smx 



sin 2 x 



xPr^ +1) (cos 



We note 

qM = ^W 2 

Then, we have 



x ) -2(£ + 3)cotxP L 1i +1) (cos X )] .(B.31) 

(B.32) 



J^ D2P l+i ( cos X)-P L +i (cosx) 





^d 2 


d 


sinx 


dx 2 


+ cotx-r- 
dx 




~d 2 


d 


sinx 


dx 2 


+ cotx-r- 



(£+1): 

sin 2 x 



^ +1) (cosx) 



+ sinx [2 - (L + 1)(L + 2)] p-^ +1) (cosx) 
« M 2 sin xPr ( ^ +1) (cos x). 



(£+1): 

sin 2 x 



(B.33) 
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Hence, we get 

Z( X ) = 2(£ + 3)D iq2 (x) + sm X Pi {£+1 \cosx)- (B.34) 
In a similar manner, in the M — > limit, we find 

M 2 W(x) - 2(£ + 3) M 2 q 2 ( X ) cot x ~ M 2 sin x-Pi~^ +1 ^(cos x)- (B.35) 
Hence, we get 

H/( X ) = 2(£ + 3)g 2 (x)cotx + sin X Pr (m) (cosx). (B.36) 



B.4 Summary 

In summary, we find 



A: 



Oto 

XX 



I- 1 
sinx 

92 (x) 
2(£ - 



<zi(x)ir, 

1 



sin 3 x 



K xi - 



k oe<r 



1 



2(^-1) /d 



£(£ + 2 
+e*P 2 ~ 
1 



^(£ + 2) 



2 



smx \dx 
(cos x) 
d_ 
dx 



+ cotx qi(x) 



+ cotx I 92 (x) 



+^3 P i ( cos x) 



No- 



where 



9i(x) = 



92 (x) 



■ 



1 



£- 1 
3a — 5 d 



29 2 (x)cotx + ^|Pr (m) (cosx) 



^-i)^^ (m) (cosx) 



to 



15 &/ 
1 



(B.37) 



15 dv 

Pr ( ' +1) (cosx) - cosxiY (m) (cosx)] } , (B.38) 
P-^Coosx) 

i/=2 

P 2 - (m) (cosx) - cosxPr (m) (cosx)l • (B.39) 



3(£ + 3) 
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